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DIGITAL MODULATION SIGNAL RECEIVER WITH ADAPTIVE CHANNEL 
EQUALIZATION EMPLOYING DISCRETE FOURIER TRANSFORMS 

This application is filed under 35 U.S.G 1 1 1(a) claiming pursuant to 35 U.S.G 
1 19(e)(1) benefit of the filing date of provisional U. S. patent application serial No, 
60/193,301 filed March 30, 2000, pursuant to 35 U.S.C. 111(b). 

The invention relates to equalization methods and apparatus used for determining 
the weighting coefficients of adaptive filtering used for equalization and echo 
cancellation in digital communications receivers for receiving digitally modulated 
single-carrier signals, as exemplified by receivers for digital television (DTV) signals as 
broadcast in the United States of America. 

BACKGROUND OF INVENTION 

Such broadcasting has been done in accordance with a Digital Television 
Standard published in 1995 by the Advanced Television Systems Committee (ATSC) as 
Document A/53. 

The component of the broadcast DTV signal to which the receiver synchronizes 
its operations is called the principal signal, and the principal signal is usually the signal 
received directly over the shortest transmission path. Therefore the multipath signal 
components of the broadcast TV signal received over other paths are usually delayed 
with respect to the principal signal and appear as lagging ghost signals. It is possible 
however, that the direct or shortest path signal is not the signal to which the receiver 
synchronizes. When the receiver synchronizes its operations to a (longer path) signal that 
is delayed respective to the direct signal, there will be a leading ghost signal caused by 
the direct signal, or there will a plurality of leading signals caused by the direct signal and 
other reflected signals of lesser delay than the signal to which the receiver synchronizes. 
In the DTV art ghost signals are customarily referred to as "echoes", because of their 
similarity to echoes in transmission lines that are terminated other than with their 
characteristic impedance. The leading ghost signals are referred to as "pre-echoes t! , and 
the lagging ghost signals are referred to as "post-echoes". The ghost signals or echoes 
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vary in number, amplitude and delay time from location to location and from channel to 
channel at a given location. Post-echoes with significant energy have been reported as 
being delayed from the reference signal by as many as sixty microseconds. Pre-echoes 
with significant energy have been reported leading the reference signal by as many as 
5 thirty microseconds. This 90-microsecond or so possible range of echoes is appreciably 
wider than television receiver designers generally supposed until midyear 2000. 

The adaptive filtering used for channel equalization and echo cancellation in 
receivers for DTV signals broadcast in accordance with the ATSC standard has generally 
been a transversal filter the kernel of which is tapped at symbol intervals for 

1 0 implementing what is termed "synchronous equalization". Synchronous equalization of 
real-only received signal is attempted, by adjusting the phase of received signal to 
minimize any imaginary component of the received signal. The attempt at adjustment is 
often made using decision-feedback methods based on the response of the adaptive 
filtering. If synchronous equalization is employed, the received signal is under-sampled 

1 5 when phase modulation of the received signal occurs during multipath reception, so 
rapidly changing phase modulation cannot be tracked by the adaptive filtering. Until 
midyear 2000, the range of echoes with appreciable energy was believed to extend from 
pre-echoes advanced by no more than three or four microseconds to post-echoes delayed 
as much as forty microseconds or so. Accordingly, the kernel of the adaptive filtering 

20 used for synchronous equalization of real-only received signal in most designs had only 
about 500 taps at symbol-epoch intervals. 

It is known generally in digital communications receiver design that "fractional 
equalization", in which the adaptive filtering kernel has taps at less-than-symbol-epoch 
intervals outperforms synchronous equalization when multipath reception conditions 

25 obtain. Page 535 of Data Communications Principles by Gitlin, Hayes & Weinstein, 
published in 1992 by Plenum Press of New York, indicates that a (3/4)-symbol-epoch 
fractional equalizer performs substantially as well as the (l/2)-symbol-epoch fractional 
equalizer known without doubt to employ adequate sampling. Fractional equalization 
costs more in die area, of course, supposing the adaptive filtering to be constructed in a 

30 monolithic integrated circuit. A (3/4)-symbol-epoch fractional equalizer has one-third 
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more taps than a synchronous equalizer, or more than 640 taps. This appears to have 
deterred fractional equalization being used in DTV signal receivers designed for use in 
homes and made available for purchase before April 2001 . 

Training signals can be used for determining the weighting coefficients of 
adaptive filtering used for channel equalization and echo suppression in receivers for 
television signals. For example, a subcommittee of the ATSC approved ghost- 
cancellation reference (GCR) signals that incorporated Bessel chirp signals on pedestals 
in the 28 th horizontal scan lines of image fields of analog television signals broadcast in 
accordance with a National Television System Committee (NTSC) standard used in the 
United States of America. The Bessel chirps in the later of the two image fields of each 
television frame are of opposite sense of polarity from the Bessel chirps in the earlier of 
the two image fields. This supports the receiver combining Bessel chirps from an 
even-numbered plurality of consecutive image fields to suppress pedestal, horizontal 
synchronizing pulse, color burst and porch information accompanying the Bessel chirps 
to generate a separated GCR signal together with the echoes, or ghosts, thereof. The 
separated GCR signal is employed as a training signal from which the weighting 
coefficients of adaptive filtering used for equalization and echo-suppression are 
determined. This determination is conveniently made by digitizing the GCR signal as 
separated from received NTSC analog television signal, determining the discrete Fourier 
transform (DFT) of successive samples of the digitized separated GCR signal and its 
ghosts as received off-the-air, and dividing that DFT term-by-corresponding term by the 
DFT of successive samples of a digitized ghost-free GCR signal as known a priori and 
stored in read-only memory at the receiver. The result of this term-by-term division is a 
DFT characterizing the actual reception channel as to its frequency response. This DFT 
is divided term-by-corresponding-term into a DFT characterizing the ideal reception 
channel as to its frequency response, to determine the DFT of the system function in the 
frequency domain of the adaptive filtering to be used for equalization and 
echo-suppression. This last DFT in the frequency domain is subjected to an inverse 
discrete Fourier transform (IDFT) procedure to generate the weighting coefficients of the 
adaptive filtering kernel in the time domain. These are closed-form computations of the 
weighting coefficients of adaptive filtering used for equalization and echo-suppression. 
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A definite solution as to the value of each weighting coefficient in the filter kernel is 
directly obtained, without open-form computations for successively approximating that 
value with reduced error over time. 

Fourier transforms in general and DFTs in particular are known to have an 
5 interesting property, which is not exploited in the adaptive filtering procedure described 
in the previous paragraph. A shift of the original data within the transform window is 
reflected solely in a change in the phasings of the transform coefficients and not in their 
amplitudes. This "Fourier transform shift theorem" was propounded by E. O. Brigham in 
"The Fast Fourier Transform" published in 1974 by Prentice-Hall of Englewood Cliffs, 
1 0 New Jersey. This property underlies fast computation of DFTs of continuous data 

streams by methods such as those described by K. B. Welles II and R. L Hartley in their 
U. S. patent No. 4,972,358 issued 20 November 1990 and titled "Computation of Discrete 
Fourier Transform Using Recursive Techniques 55 . 

Blind-equalization methods have been resorted to for determining the weighting 

1 5 coefficients of adaptive filtering used for equalization and echo-suppression in receivers 
for ATSC digital television (DT V) signals, because the ATSC DTV signal does not 
provide a good training signal for procedures similar to those described in the preceding 
paragraph. The data field synchronizing (DFS) signals specified by ATSC Document 
A/53 each include a PN51 1 pseudo-random noise signal and a triple PN63 

20 pseudo-random noise signal. At the time Document A/53 was published, these 

pseudo-random noise (PN) signals were envisioned as being used as training signals for 
adaptive channel-equalization and echo-suppression filtering. The design of the DFS 
signal does not avoid the PN signals being overlapped by the echoes of previous data in 
the DTV signal that have significant energy, however, nor does the design avoid some 

25 echoes of these PN signals that have significant energy overlapping subsequent data in 
the DTV signal Consequently, even though the PN signals have auto-correlation 
functions that might suit them for match filtering, the DFS signals have not proven in 
actual practice to be satisfactory training signals for adaptive equalization and 
echo-suppression filtering, because echoes of the PN signals are not readily 

30 distinguishable from other data and their echoes. 
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Therefore, data-directed methods have been resorted to for computing on a 
continuing basis the weighting coefficients of the adaptive filtering used for equalization 
and echo-suppression in receivers for ATSC DTV signals. The approaches usually are 
auto-regressive spectral analyses, which are generally described as follows. The actual 
5 response of the adaptive filtering to received signal is data-sliced, or quantized, to 
generate an estimate of the symbols actually transmitted. The actual response of the 
adaptive filtering to received signal is compared symbol-epoch-by-symbol-epoch with 
the estimates of the symbols actually transmitted to generate an error signal to be used in 
a decision-feedback procedure for calculating the weighting coefficients of the adaptive 
1 0 filter. The decision-feedback procedure uses one of a variety of known algorithms that 
operate on a successive approximation basis. 

LMS-gradient algorithms used in data-directed equalization methods are quite 
slow and are prone to stalling at local minima in the decision-feedback error signal, rather 
than continuing to the minimal decision-feedback error signal condition overall. 
1 5 Initialization of the adaptive filtering used for equalization and ghost cancellation 
normally takes a second or so after the reception channel is changed, which makes 
channel surfing difficult. Some receivers store weighting coefficients from the last time a 
channel was tuned to, to furnish a starting point for initializing the adaptive filtering. 

LMS-gradient algorithms adapt slowly, so rapidly occurring changes in multipath 
20 conditions cannot be followed. This will at times lead to the weighting coefficients being 
completely erroneous for changed multipath condition, causing data slicing errors 
frequently enough that the error-correction capabilities of the system are overwhelmed. 
In some instances the multipath conditions do not return to a previous state which the 
weighting conditions are reasonably correct for, to cause data slicing errors to be 
25 infrequent enough that they can be corrected. Then, there is a second or so interval after 
rapid changes in multipath conditions pass before the LMS-gradient algorithm can 
re-initialize the weighting coefficients of the adaptive filtering. Continuously changing 
multipath conditions can cause a loss of tracking in which data slicing errors too frequent 
to be corrected persist over protracted intervals many seconds long. 
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Alternatively, recursive least squares (RLS) filter adaptation methods are the 
data-directed equalization methods used for computing the weighting coefficients the 
adaptive filtering a DTV receiver uses for equalization and ghost-cancellation. If 
signal-to-noise conditions are high, the RLS algorithm converges in about 2M+2 
5 iterations for small error signals, where M+l is the number of taps in the kernel of the 
adaptive filter, which is typically about an order of (binary) magnitude faster than 
LMS-gradient algorithms converge, Such faster convergence would appear to reduce the 
chances for loss of adaptive filter tracking owing to dynamic multipath conditions and for 
such loss being of protracted duration. However, the tracking performance is influenced 

10 not only by the rate of convergence (which is a transient characteristic) but also by 
fluctuation of the steady-state performance of the algorithm as influenced by 
measurement and algorithm noise. With both algorithms tuned to minimize the 
misadjustment of the filter response by a proper optimization of their forgetting rates, the 
LMS algorithm exhibits tracking performance superior to that of the RLS algorithm. 

1 5 Moreover, dynamic multipath conditions tend to be more troublesome during 

weak-signal reception where the convergence of the RLS algorithm is not so much 
greater than that of the LMS-gradient algorithms. 

RLS methods generally involve a computational cost that increases about as the 
square of the number of taps contained in the adaptive filter. This is a prohibitively high 

20 cost for a DTV receiver designed for use in homes. The fast transversal filters (FTF) 

algorithms realize the RLS solution with a computational cost that increases only linearly 
with the number of taps contained in the adaptive filter, as in the LMS-gradient 
algorithms. The computational cost of the FTF algorithms is at least four times larger 
than that for the LMS-gradient algorithms, however, with division calculations being 

25 required For an adaptive filter having an (M+l)-tap kernel, the LMS-gradient algorithms 
require about 2M+1 multiplications and 2M additions/subtractions^ with no divisions 
being required. The FTF algorithms require at least 7M+12 multiplications, 4 divisions, 
and 6M+3 additions/subtractions, with additional computation being required to avoid 
long-term instabilities. 
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Like the LMS-gradient algorithms, the FTF algorithms can suffer from stalling at 
local minima in the error signal, but methods are known for preventing such stalling. The 
more intractable problem with FTF algorithms is a potential explosive instability arising 
from word-length limitations in the weighting-coefficient computer causing accumulated 
5 rounding errors. One method that has been used for avoiding this explosive instability is 
to evaluate error signals developed by comparing the results of alternative ways of 
calculating the FTF algorithm, the differences in results being attributed to errors 
introduced by rounding off to accommodate word-length limitations in the 
weighting-coefficient computer. Such methods increase computational complexity by 

10 another 1 5% or so. Another method that has been used for avoiding this explosive 

instability is periodically restarting the FTF algorithm, with an LMS-gradient algorithm 
taking over in the interim proceeding from the weighting coefficients the FTF algorithm 
has computed. The LMS-gradient algorithm eliminates accumulated errors in the 
coefficients, so the FTF algorithm can resume its calculations without accumulated 

1 5 round-off errors. The hand-off to the LMS-gradient algorithm also increases 
computational complexity. 

The adaptive filters that have previously been used for equalization and 
echo-suppression in receivers for ATSC DTV signals are tracking filters that perform 
auto-regressive spectral analyses. When deep fading occurs suddenly during dynamic 

20 muitipath reception conditions, there is a tendency for tracking to be lost. Tracking of the 
adaptive filtering may be impossible to recover unless a complete re-initialization of its 
weighting coefficients can be accomplished before there is another sudden change in 
dynamic muitipath reception conditions. DTV receiver designers have attempted to solve 
the loss of tracking problem by improving the tracking rate of the adaptive filtering. This 

25 attempts to prevent loss of tracking, rather than dealing with the problem of re-initializing 
filter coefficients "instantly" when tracking is lost. The problem is that objects moving at 
fairly low velocity can at times interrupt reception from one of two reception paths of 
similar strength to cause tremendous changes in signal phase that are nearly 
"instantaneous". Accordingly there will always be times, hopefully rare, when adaptive 

30 tracking of the adaptive filter will be lost. The question then circles back to how rapidly 
the adaptive filter coefficients can be re-initialized after tracking is lost, particularly when 
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reception conditions do not change back to a previous state. A few seconds is simply too 
long, particularly since audio as well as video is lost. 

Auto-regressive spectral analyses are handicapped in regard to how rapidly the 
weighting coefficients of the adaptive filtering can be re-initialized after tracking is lost 
5 The ultimate limit on how rapidly re-initialization is possible after loss of tracking is 
determined by the desideratum that re-initialization of the weighting coefficients should 
be deferred until the data on which their computation is based no longer appreciably 
affect the adaptive filtering response. Such deferral is necessary in order to assure 
stability of the feedback loops in which error signal is derived from the adaptive filtering 
1 0 response to support computation of the weighting coefficients. 

The speed of re-initialization is limited not only by the latency of the adaptive 
filtering used for equalization and echo-suppression, however, but additionally by the 
latency associated with computation of the weighting coefficients from the 
decision-feedback error signals. The updates of weighting coefficients are not computed 

1 5 in parallel in the auto-regressive methods, but are computed seriatim based on 

minimizing the decision-feedback error signals over time. The updates of the weighting 
coefficients can be applied to the adaptive filter serially as they are computed, but most 
designs avoid undesirable reverse-time effects by periodically applying the updates of the 
weighting coefficients to the adaptive filter so as to update the weighting coefficients in 

20 the entire kernel simultaneously in a technique known as "block updating". The serial, 
rather than parallel, computation of updates for the weighting coefficients slows the 
adaptation of the filtering used for equalization and echo-suppression, particularly when 
block updating of the kernel weighting coefficients is done. 

Furthermore, the auto-regressive methods employ open-form computation for 
25 continually adjusting the kernel weighting coefficients of the adaptive filtering using 

successive approximation techniques, rather than computing those weighting coefficients 
outright using closed-form computation. These open-form computations convolve the 
adaptive filter kernel with several blocks of signal samples, rather than a single block of 
samples. This tends to make initial computation of an (M+l)-tap kernel an operating 
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method with steps some (3M+2) times N in number. This sort of operating method tends 
to take considerable time to perform, since these steps N(3M+2) in number have to be 
performed at the normal sampling rate through the adaptive filter 

The factor (3M+2) arises from (2M+1) sample epochs being required for 
5 convolving (M+ 1) samples of input signal to the adaptive filter with the (M4-l)-tap kernel 
to generate the (M+l) samples of decision-feedback signal required for updating all the 
weighting coefficients in the kernel, and from (2M+1) sample epochs being required for 
convolving the (M+l) samples of the decision-feedback signal gradient with the 
(M-U)-tap kernel to generate the updates of ail the weighting coefficients in the kernel. 
10 The latter convolution procedure is presumed to commence next sample epoch after the 
first sample of decision-feedback signal is generated by the former convolution 
procedure, so the latter convolution procedure overlaps the former convolution procedure 
over M sample epochs. 

The factor N is reciprocally related to an attenuation factor that is introduced into 

1 5 each successive one of the computations to update the weighting coefficients applied to a 
respective one of the kernel taps. The factor N is introduced into these computations so 
that the optimum value of each weighting coefficient is approached or reached through 
successive approximation. That is, the weighting coefficients are generated in the 
auto-regressive adaptive-filtering methods by long-term accumulations of 

20 decision-feedback error signal energy supplied in small increments and small decrements. 
These procedures result in lowpass recursive filtering of each weighting coefficient in the 
adaptive filtering kernel, which suppresses the effects of noise on the computation of the 
weighting coefficients. The noise comprises quantization noise generated in the digital 
portion of the receiver, Johnson noise from the analog portion of the receiver, impulse 

25 noise in the reception channel and possibly co-channel interference from NTSC analog 
television signals or from other ATSC DTV signals. This form of lowpass filtering of 
each weighting coefficient in the adaptive filtering kernel to suppress the effects of noise 
on its computation, which lowpass filtering is an integral part of the accumulation 
procedure that implements the computation, has a severe shortcoming as compared to 

30 lowpass filtering of each weighting coefficient and its updates after their generation. 
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That is, since the weighting coefficients are computed seriatim in the auto-regressive 
adaptive-filtering methods, the number N of computations of each weighting coefficient 
that are averaged to arrive at the final value of that weighting coefficient actually 
employed in the adaptive filtering kernel appears as a multiplicative factor in determining 
5 the number of sample epochs required for completing initialization of the kernel 

weighting coefficients. If lowpass filtering of each weighting coefficient and its updates 
after their generation were used instead, the number N of computations of each weighting 
coefficient that are averaged to arrive at the final value of that weighting coefficient 
actually employed in the adaptive filtering kernel would simply add N sample epochs to 
10 the number of sample epochs required for completing initialization of the kernel 
weighting coefficients. 

This specification discloses a novel equalization method for adapting the kernel 
weighting coefficients of the filtering used for equalization and echo-suppression in DTV 
signal receivers. This novel equalization method uses closed-form computations to 

1 5 obtain directly a complete update of the value of each weighting coefficient in the filter 
kernel, without open-form computations for successively approximating that value with 
reduced error over time. Closed-form computations that perform convolution of 
(M+l)-sample terms using discrete-Fourier-transform (DFT) procedures and convert the 
results obtained in the frequency domain back to the time domain using 

20 inverse-discrete-Fourier-transform (I-DFT) procedures can be performed in as few as 

2(M+1) log 2 (M+l) sample epochs. The computations to determine the complete update 
of each weighting coefficient in the filter kernel are performed in parallel After the 
complete updates are computed in parallel, all the weighting coefficients in the kernel of 
the adaptive filtering used for equalization and echo-suppression are simultaneously 

25 updated in full 

While the novel equalization method includes steps of lowpass filtering each 
weighting coefficient included in the adaptive filtering kernel, in order to suppress the 
effects of noise on the computation of the weighting coefficients, the lowpass filtering is 
performed after full updates of coefficient values have been computed So ? lowpass 
30 filtering of N successive values of each of the parallelly computed weighting coefficients 
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adds only N sample epochs to the number of sample epochs required for completing 
initialization of the kernel weighting coefficients. N does not appear as a multiplicative 
factor in determining the number of sample epochs required for completing initialization 
of the kernel weighting coefficients. 

5 Complete updating of the kernel weighting coefficients by the novel 

equalization method can initially be done in N + 2(M+1) log 2 (M+l) sample epochs. 
During initialization or re-initialization N can be made to be zero, or to be otherwise 
smaller than it is during continuing operation thereafter. 

N is usually an integral power of two, better to implement lowpass filtering that 
10 employs a tree of digital adders. Typical values are 32, 64, 128 or 256, with higher 
values of N being favored in auto-regressive adaptive filtering methods. U. S. patent 
application serial No. 60/193,301 indicates that (M+l) would preferably have a value at 
least 5 12 for representing a time period of 30-50 microseconds duration. That patent 
application further indicates that value of (M+l) = 1024 would be preferred for a 
15 fractional equalizer using twice-baud-rate sampling and having its weighting coefficients 
adapted by the novel equalization method described in this specification. IX S. patent 
application serial No. 60/193,301 indicates that assured initialization or 
re-initialization takes about two milliseconds with preferred forms of the novel 
equalization method described in its specification, as compared to the 
20 25-millisecond fastest initialization or re-initialization time claimed for auto-regressive 
blind-equalization methods. The novel equalization method described in U. S. patent 
application serial No. 60/193,301 and in this specification is never subject to the stalling 
problems that can afflict auto-regressive blind-equalization methods. 

In order for DTV receiver designs employing a synchronous equalizer to 
25 accommodate the 90-microsecond-duratkm echo range that was publicly disclosed in 
2000, (M+l) preferably has a value 2 P where P is eleven or more. In order that a DTV 
receiver design employing a fractional equalizer using twice-baud-rate sampling can 
accommodate the 90-microsecond-duration echo range, (M+l) preferably has a value 
2 (P+l) . Presuming that the duration of the window for the DFT is increased from a 47.6 
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microsecond duration to 190,3 microseconds, assured initialization or re-initialization 
will take another millisecond or so longer. 

SUMMARY OF THE INVENTION 

The invention generally concerns novel equalization techniques for calculating 
5 the system characteristic of the adaptive filtering used for equalization and 

echo-suppression in a digital communications receiver, such as one used for receiving 
over-the-air broadcast digital television signals, and the receiver circuitry for 
implementing these novel equalization techniques. In these novel equalization 
techniques, the system characteristic of the adaptive filtering is calculated from the 
1 0 discrete Fourier transform of the input signal supplied to the adaptive filtering and from 
the discrete Fourier transform of the transmitted signal as estimated in the receiver, 

A species of this novel equalization technique computes the weighting 
coefficients for adaptive transversal filtering used for equalization and echo-suppression 
of ATSC DTV signals in receivers constructed in accordance with the inventions 

15 described in this specification. In this technique the computation of the weighting 

coefficients of the filtering is based on channel characterization that is computed in the 
following manner, assuming the adaptive filter receives baseband ATSC DTV signal as 
input signal. The adaptive filtering response is supplied to a symbol decoder, such as a 
data sheer or a Viterbi trellis decoder, which estimates the symbols transmitted from the 

20 transmitter. The discrete Fourier transform (DFT) of a large plurality P in number of 
successive samples of the transmitted symbols as so estimated is computed A 
corresponding large plurality P in number of successive samples of the signal supplied to 
the adaptive filtering for equalization and echo-suppression is delayed, to compensate for 
the latency of the adaptive filtering and the symbol decoder cascaded thereafter, and the 

25 DFT of the delayed plurality of successive samples is calculated. The DFT of the 

plurality P in number of successive samples of the transmitted symbols, as estimated, is 
multiplied term-by-corresponding term by a DFT characterizing the ideal reception 
channel with Nyquist roll-off for minimizing inter-symbol interference (ISI), which latter 
DFT is read from read-only memory in the receiver. The DFT defined by the products 
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resulting from this term-by-corresponding term multiplication describes the DFT of the 
desired adaptive filter response to the adaptive filtering input signal. The DFT of the 
adaptive filtering response to a unit impulse input signal that is needed for equalization 
and echo-suppression of the actual reception channel response is computed by dividing 
5 each term of the DFT of the desired adaptive filter response to a plurality of successive 
samples of adaptive filtering input signal by the corresponding term of the DFT of the 
plurality of successive samples of adaptive filtering input signal, as delayed to 
compensate for the latent delay of the adaptive filter. The inverse discrete Fourier 
transform (IDFT) of the DFT of the adaptive filtering response specifies in the time 
10 domain the kernel of the adaptive transversal filtering for equalization and 

echo-suppression. Of course, owing to the commutative and associative properties of 
multiplication and division processes, at least one variation of this procedure exists and 
produces an equivalent result 

Another species of the novel equalization technique is implemented with adaptive 
1 5 filtering that uses a parallel-bank analysis filter for determining the DFT of the adaptive 
filtering input signal; combining circuitry for multiplying the DFT of that input signal 
term-by-corresponding-term with the DFT of the adaptive filtering response that is 
needed for equalization and echo-suppression of the actual reception channel response, 
thereby to generate the DFT of the adaptive filtering output signal; and a synthesis filter 
20 for generating the adaptive filtering output signal from its DFT, by inverse-DFT 

computation. The adaptive filtering output signal is supplied to the symbol decoder, 
which estimates the symbols transmitted from the transmitter The DFT of these 
estimates is calculated and is then processed together with the DFT of the adaptive 
filtering input signal for generating the DFT of the adaptive filtering response that is 
25 needed for equalization and echo-suppression of the actual reception channel response. 

Passband equalization is done in accordance with further aspects of the invention, 
using procedures similar to the procedures for baseband equalization described in the 
foregoing two paragraphs. 

DESCRIPTION OF THE DRAWING 
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FIGURE 1 is a schematic diagram of a novel digital television receiver, which 
uses adaptive filtering of real-only signal as sampled at higher than baud rate, and which 
uses data slicing of the adaptive filtering response to generate estimates of actually 
transmitted VSB DTV signal symbols. 

5 FIGURE 2 is a flow chart of a novel method that can be used in the FIGURE 1 or 

FIGURE 3 DTV receiver for computing the kernel weighting coefficients of the adaptive 
filtering used for equalization and echo-suppression. 

FIGURE 3 is a schematic diagram of a novel digital television receiver, which 
uses adaptive filtering of real-only signal as sampled at higher than baud rate, and which 
10 uses Viterbi decoding of the adaptive filtering response to generate estimates of actually 
transmitted VSB DTV signal symbols. 

FIGURE 4 is a schematic diagram of a novel digital television receiver that is a 
modification of the FIGURE 1 DTV receiver. 

FIGURE 5 is a flow chart of a novel method that can be used in the FIGURE 4 or 
1 5 FIGURE 6 DTV receiver for computing the kernel weighting coefficients of the adaptive 
filtering used for equalization and echo-suppression. 

FIGURE 6 is a schematic diagram of a novel digital television receiver that is a 
modification of the FIGURE 3 DTV receiver. 

FIGURE 7 is a schematic diagram showing in detail demodulator and 
20 over-sampling analog-to-digital conversion circuitry that can be incorporated in any of 
the digital television receivers of FIGURES 1 ? 3, 4 and 6. 

FIGURE 8 is a schematic diagram of portions of a novel digital television 
receiver that employs an adaptive filter for passband equalization and echo-suppression, 
which adaptive filter has its kernel weighting coefficients computed in accordance with 
25 the invention. 

FIGURE 9 is a schematic diagram of a novel digital television receiver that is a 
modification of the FIGURE 8 DTV receiver. 
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FIGURE 10 is a schematic diagram of portions of a novel digital television 
receiver that employs an adaptive filter for complex baseband equalization and 
echo-suppression, which adaptive filter has its kernel weighting coefficients computed in 
accordance with the invention. 

5 FIGURE 1 1 is a schematic diagram of portions of another novel digital television 

receiver that employs an adaptive filter for complex baseband equalization and 
echo-suppression, which adaptive filter has its kernel weighting coefficients computed in 
accordance with the invention. 

FIGURE 12 is a schematic diagram of portions of a novel digital television 
10 receiver that employs an adaptive filter for complex passband equalization and 

echo-suppression, which adaptive filter has its kernel weighting coefficients computed in 
accordance with the invention. 

FIGURE 13 is a schematic diagram of portions of another novel digital television 
receiver that employs an adaptive filter for complex passband equalization and 
1 5 echo-suppression, which adaptive filter has its kernel weighting coefficients computed in 
accordance with the invention. 

FIGURE 14 is a schematic diagram of a modification that is made in the adaptive 
filtering in the FIGURE 1, 3, 4 and 6 DTV signal receivers in other embodiments of the 
invention. 

20 FIGURE 1 5 is a schematic diagram of a modification that is made in the adaptive 

filtering in the FIGURE 8 and 9 DTV signal receivers in other embodiments of the 
invention. 

FIGURE 16 is a schematic diagram showing how an alternative adaptive filter 
structure is employed in a modified FIGURE 1 DTV signal receiver that embodies the 
25 invention. 
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FIGURE 17 is a schematic diagram showing how an alternative adaptive filter 
structure is employed in a modified FIGURE 3 DTV signal receiver that embodies the 
invention. 

FIGURE 18 is a schematic diagram showing how an alternative adaptive filter 
5 structure is employed in a modified FIGURE 4 DTV signal receiver that embodies the 
invention. 

FIGURE 19 is a schematic diagram showing how an alternative adaptive filter 
structure is employed in a modified FIGURE 6 DTV signal receiver that embodies the 
invention. 

10 FIGURE 20 is a schematic diagram showing how an alternative adaptive filter 

structure is employed in a modified FIGURE 8 DTV signal receiver that embodies the 
invention. 

FIGURE 21 is a schematic diagram showing how an alternative adaptive filter 
structure is employed in a modified FIGURE 9 DTV signal receiver that embodies the 
15 invention. 

DETAILED DESCRIPTION 

In the FIGURE 1 digital television receiver apparatus a source 10 of 
radio-frequency vestigial-sideband digital television signal, which typically is an antenna 
for receiving that R-F VSB DTV signal as broadcast over the air from a terrestrial 

20 transmitter, supplies that R-F VSB DTV signal to a DTV receiver front-end 11. The 
DTV receiver front-end 11 typically comprises a radio-frequency (R-F) amplifier, an 
up-converter for converting the amplified R-F VSB DTV signal from the R-F amplifier to 
an ultra-high-frequency (UHF) VSB DTV signal with a carrier frequency around 920 
MHz ? a UHF intermediate-frequency amplifier, a down-converter for converting the 

25 amplified UHF VSB DTV signal from the UHF I-F amplifier to a very-high-frequency 
(VHF) VSB DTV signal in a 41-47 MHz frequency range, and an 
automatic-gain-controlled VHF intemiediate-frequency amplifier. Several alternative 
designs for the DTV receiver front-end 11 are known. In a plural-conversion receiver the 



16 



Atty. Dkt ALRL03A 



PATENT Application of J. D. McDonald et al. 



final-intermediate-frequency VSB DTV signal may repose in a frequency range lower 
than the 41-47 MHz frequency range, such as 21-27 MHz, to facilitate digitization of the 
final-intermediate-frequency VSB DTV signal followed by demodulation in the digital 
regime. Triple-conversion receivers in which down-conversion to a VHF VSB DTV 
5 signal in the 4 1 -47 MHz frequency range is followed by a further down conversion to a 
low-band intermediate-frequency band offset from zero-frequency by a megahertz or so 
are favored by some DTV receiver designers. 

FIGURE 1 shows the DTV receiver front-end 11 connected for supplying 
final-I-F VSB DTV signal to demodulator and over-sampling analog-to-digital 

10 conversion circuitry 12. The circuitry 12 can take any one of several known forms. The 
final-I-F VSB DTV signal can be digitized prior to demodulation in the digital regime, as 
noted in the previous paragraph. Alternatively, the final-I-F VSB DTV signal is not 
digitized and is instead demodulated in the analog regime, after which the baseband 
demodulation results are digitized. The order in which the demodulation and 

1 5 analog-to-digital conversion of the VSB DTV signal is done in the circuitry 12 is not of 
material consequence to the FIGURE 1 embodiment of the invention, in which only the 
real component of demodulated signal is supplied to an adaptive filter 13 for equalization 
and echo-suppression. The digital samples of the real component of demodulated signal 
supplied to the adaptive filter 13 should be generated at a sampling rate higher than baud 

20 rate, or symbol rate, of the transmitted DTV signal. Otherwise, the VSB DTV signal as 
received with dynamic multipath distortion will be under-sampled, because of the angular 
modulation caused by dynamic multipath distortion. Generating the digital samples of 
the real component of demodulated signal supplied to the adaptive filter 13 at twice baud 
rate will assure that the Nyquist criterion for adequate sampling is met irrespective of the 

25 changes in the phasing of the received DTV signal owing to multipath distortion. 

The I-F amplifiers of the FIGURE 1 DTV receiver limit the amount of 
excess-bandwidth available for phase-modulation sidebands, so sampling at rates 
somewhat lower than twice baud rate is probably sufficient The experience of others in 
digital communications receivers indicates that sampling at four-thirds baud rate to 
30 implement (3/4>symbol-epoch fractional equalization provides substantially as good 
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performance as sampling at twice baud rate to implement (l/2>symbol~epoch fractional 
equalization. See page 535 of Data Communications Principles by Gillin, Hayes & 
Weinstein, copyright 1992 to Plenum Press of New York and London. 

FIGURE 1 shows the response of the adaptive filter 13 used for equalization and 

5 echo-suppression being supplied to decimation fitter circuitry 14 for resampling to baud 
rate. FIGURE 1 shows the decimation filter circuitry 14 connected for supplying its 
response to a trellis decoder 15 of the twelve-phase type conventionally employed in an 
ATSC VSB DTV receiver. The trellis decoder 15 supplies data resulting from symbol 
decoding to a data de-interleaver 16, which undoes the convolutionai interleaving done at 

10 the ATSC VSB DTV over-the-air (OTA) broadcast transmitter, so as to disperse 

extended-duration noise bursts that result from excitation of the I-F amplifiers by impulse 
noise. The data de-interleaver 16 is connected to supply de-interleaved data in bytes to 
circuitry 17 for correcting errors in response to Reed-Solomon forward error correction 
coding contained in the de-interleaved data. The circuitry 17 is connected to supply 

1 5 error-corrected data to a data de-randomizer 18, which regenerates the stream of data 
packets supplied to the ATSC VSB DTV broadcast transmitter for over-the-air 
broadcasting. This stream of data packets is supplied to the rest 19 of the FIGURE 1 
ATSC VSB DTV receiver, which selects data packets descriptive of compressed video to 
an MPEG-2 decoder, selects data packets descriptive of audio to an AC-3 decoder, and 

20 performs the other functions of processing the video and audio for presentation or for 
recording to be presented at a later time. In sum, elements 15 - 19 function similarly to 
the way such elements function in other ATSC VSB DTV receivers previously known. 

Consider now how the weighting coefficients in the kernel of the adaptive filter 
13 are determined in accordance with the invention. The decimation filter circuitry 14 

25 response is supplied to a data sheer 20, which detects data slices by a "hard"-decision 
symbol decoding method. The data sheer 20 supplies its "hard"-decision symbol 
decoding results as input addressing to a read-only memory 21 that recodes these 
decoding results using ideal symbol code levels. The symbol coder ROM 21 output 
signal provides estimates of the symbols that were actually transmitted by the VSB DTV 

30 OTA broadcast transmitter. The boundaries of the data slicing by the data slicer 20 
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preferably take into account the pilot carrier accompanying the received symbols. That 
is, the boundaries of the data slicing by the data slicer 20 are asymmetrically disposed 
relative to the zero-energy axis for detecting normalized modulation levels of -5.75, 
-3.75, -1.75, +0.25, +2.25, +4.25, +6.25, + 8.25 nominal values. Such data slicing 
5 facilitates arrangements for extracting automatic frequency and phase control (AFPC) 
information from the adaptive filter 13 response for controlling the local oscillations used 
to synchrodyne the VSB DTV signal to base band in synchronous demodulator portions 
of the circuitry 12. 

However, in other embodiments of the invention the low-frequency portion of the 
10 adaptive filter 13 input signal near zero frequency may be removed in the demodulator 
and ADC circuitry 12. In such case the boundaries of the data slicing by the data slicer 
20 take into account the received symbols not being accompanied by pilot carrier and are 
symmetrically disposed relative to the zero-energy axis for detecting normalized 
modulation levels of -7, -5, -3, -1, +1, +3, +5, +7 nominal values. 

15 The estimates of actually transmitted symbols provided by the symbol coder 

ROM 21 are supplied to interpolation filter circuitry 22 for re-sampling back to the 
sampling rate used in the adaptive filter 13. DFT computation circuitry 23 computes the 
discrete Fourier transform (DFT) of the estimation of the signal that the VSB DTV OTA 
broadcast transmitter actually transmitted to the FIGURE 1 VSB DTV receiver, as 

20 supplied by the interpolation filter circuitry 22. 

The estimation of the signal that the VSB DTV OTA broadcast transmitter 
actually transmitted to the FIGURE 1 VSB DTV receiver, as supplied by the 
interpolation filter circuitry 22, is delayed respective to the actually received signal 
supplied from the demodulator and over-sampling analog-to-digital conversion circuitry 
25 12 to the adaptive filter 13 as input signal thereto. The input signal to the adaptive filter 
13 is subjected to compensatory digital delay 24 equal to the combined latent delays of 
the adaptive filter 13, the decimation filter circuitry 14, the data slicer 20, the symbol 
coder ROM 21 and the interpolation filter circuitry 22. DFT computation circuitry 25 
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computes the discrete Fourier transform (DFT) of the delayed input signal to the adaptive 
filter 13 that the compensatory digital delay 24 provides. 

Combining circuitry 26 combines the DFTs generated by the DFT computation 
circuitry 23 and the DFT computation circuitry 25 with the DFT of ideal channel 
5 response as supplied from read-only memory 27 to generate the DFT of the adaptive filter 
13 kernel weighting coefficients. This ideal channel response is a lowpass filter response 
with Nyquist slope raised-cosine roll-off to minimize intersymbol interference (ISI), This 
lowpass filter response has flat low-frequency response near zero frequency in some 
embodiments of the invention. 

10 In other embodiments of the invention this lowpass filter response has boosted 

low-frequency response near zero frequency, This latter arrangement boosts the 
low-frequency response of the adaptive filter 13 near zero frequency as well, which 
facilitates arrangements for extracting automatic frequency and phase control (AFPC) 
information from the adaptive filter 13 response for controlling the local oscillations used 

15 to synchrodyne the VSB DTV signal to base band in synchronous demodulator portions 
of the circuitry 12. The low-frequency boost of the adaptive filter 13 response near zero 
frequency is then removed by suitable finite-impulse-response (FIR) digital filtering to 
generate a response suitable for application as input signal to the trellis decoder 15 and to 
the data slicer 20. The latency of this FIR digital filter as well as the combined latent 

20 delays of the adaptive filter 13, the decimation filter circuitry 14, the data slicer 20, the 
symbol coder ROM 21 and the interpolation filter circuitry 22 is then compensated for by 
the compensatory digital delay 24. 

The combining circuitry 26 supplies the DFT of the adaptive filter 13 kernel 
weighting coefficients that it generates to an inverse discrete Fourier transform computer 
25 28. The I-DFT computer 28 generates the time-domain adaptive filter 13 kernel 

weighting coefficients corresponding to the DFT supplied by the combining circuitry 26. 
These adaptive filter 13 kernel weighting coefficients in the time domain are supplied to 
respective ones of a bank 29 of digital lowpass filters. In some designs the digital 
lowpass filters are of finite-impulse-response (FIR) type; in other designs the digital 
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lowpass filters are of infinite-impulse-response (IIR) type. This bank 29 of digital 
lowpass filters suppresses fluctuations in the values of the adaptive filter 13 kernel 
weighting coefficients caused by noise. The noise comprises quantization noise 
generated in the digital portion of the receiver, Johnson noise from the analog portion of 
5 the receiver, impulse noise in the reception channel and possibly co-channel interference 
from NTSC analog television signals or from other ATSC DTV signals, The bank 29 of 
digital lowpass filters supplies the lowpass-filtered adaptive filter 13 kernel weighting 
coefficients to coefficient registers in the adaptive filter 13 to be used as multiplier input 
signals by digital multipliers in the adaptive filter 13. The adaptive filter 13 is presumed 
10 to have the canonic transversal filter structure in which a tapped digital delay line 
supplies differently delayed tap responses to the adaptive filter 13 input signal for 
application to the digital multipliers as multiplicand input signals, and in which the 
products from the digital multipliers are summed to generate the adaptive filter 13 output 
signal. 

1 5 FIGURE 2 is a flow chart showing in more detail a method 30 comprising steps 

31 - 38 by which the adaptive filter 13 kernel weighting coefficients can be calculated in 
the FIGURE 1 DTV receiver. Initially, the DFT of the delayed adaptive-filter 13 input 
signal supplied from the compensatory digital delay 24 is computed by the DFT 
computer 25 in the step 31 using a fast Fourier transform algorithm of the well-known 

20 Cooley-Tukey type, and the estimation of the actually transmitted signal as re-sampled by 
the interpolation filter 22 is computed by the DFT computer 23 in the step 32. Following 
the step 31, the DFT of the delayed adaptive-filter 13 input signal is continuously 
re-computed by the DFT computer 25 in a repeated updating step 33 using a Fourier-shift 
method such as disclosed by K. B. Welles II and R. I. Hartley in their U. S. patent No. 

25 4,972,358, Following the step 32 using a fast Fourier transform algorithm of the 

well-known Cooley-Tukey type for initial computation of the DFT of the interpolation 
filter 22 response, the DFT of the interpolation filter 22 response is continuously 
re-computed in a repeated updating step 34 performed by the DFT computer 23 using the 
same Fourier-shift method that the DFT computer 25 uses in performing the repeated 

30 updating step 33. 
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The combiner 26 performs the steps 35 and 36 of the method 30. In the step 35, 
the DFT of the interpolation filter 22 response as computed by the steps 32 and 34 is 
multiplied term-by-corresponding-term by the DFT of the frequency spectrum of the 
system function for the ideal reception channel, which latter DFT is extracted from the 

5 ROM 27. The step 35 generates a DFT of Nyquist filtering response to the estimation of 
the baseband modulating signal sent from the transmitter. In the step 36, the DFT 
generated in the step 35 is divided term-by-corresponding-term by the DFT of the 
delayed adaptive-filter 13 input signal, as computed by the steps 31 and 33. The step 36 
generates a DFT that defines in the frequency domain the ideal adaptive filter 13 kernel 

10 weighting coefficients in the time domain. Divisors that are so small they would cause 
out-of-range quotients should be augmented in the step 36 division procedures. The 
multiplications in step 35 and the divisions in step 36 can be performed logarithmically 
using tables stored in read-only memory, for example. In an alternative design approach 
the ROM 27 is dispensed with, and the steps 35 and 36 are merged into a single step 

1 5 implemented by read-only memory. 

In step 37 the computer 28 computes the inverse-DFT of this DFT of the ideal 
adaptive filter 13 kernel weighting coefficients in the frequency domain to obtain the 
weighting coefficients for the adaptive filter 13 kernel. During initialization or 
re-initialization, these weighting coefficients can be directly applied to the weighting 
20 coefficients register(s) of the adaptive filter 13 to speed up initialization or 

re-initialization as much as possible. Generally, however, in the step 38 each of the 
weighting coefficients is lowpass filtered for reducing the effects of noise before being 
supplied to the weighting coefficients registers) of the adaptive filter 13, this lowpass 
filtering being done by the bank 29 of digital lowpass filters. 

25 The DFT computations using a Fourier-shift method are performed on a 

"sliding-window" in which the window for the DFT calculation advances one sample 
time each time DFTs are re-calculated. Since the March 30, 2000, filing of provisional 
U. S. patent application serial No. 60/193,301, it was determined from simulations 
conducted by one of the inventors that the computation of DFTs on a sliding-window 

30 basis is not necessary, at least for reception with stationary receiving antennas. Instead, 



22 



Atty. Dkt. ALRL03A 



PATENT Application of J. D. McDonald et al. 



successive independent computations of DFTs can be done on a consecutive 
block-after-block basis. Excellent echo-suppression has been obtained using DFTs 
extracted from successive 2048-symbol-epoch blocks of baseband DTV signal and 
successive 2048-symbol-epoch blocks of estimated DTV modulating signal. The blocks 
of baseband DTV signal and the blocks of estimated DTV modulating signal had to 
contain substantially fewer symbol epochs for the equalization procedure to be started up 
from the data field synchronizing signal specified in ATSC Document A/53. If data 
fields are modified to include three data segments of a repetitive-PN1023 sequence with 
baud-rate symbols, the equalization procedure will start up using successive 
2048-symbol-epoch blocks of baseband DTV signal and successive 2048-symbol-epoch 
blocks of estimated DTV modulating signal. 

FIGURE 3 diagrams a DTV receiver that differs from the FIGURE 1 DTV 
receiver in that it does not use data slicing of the adaptive filtering response to generate, 
in accordance with the so-called "hard"-decision technique, the input addressing of the 
symbol coder ROM 21 used for estimating actually transmitted VSB DTV signal 
symbols. Instead, the FIGURE 2 DTV receiver employs Viterbi decoding of the adaptive 
filtering response response to generate, in accordance with the so-called "soff'-decision 
technique, the input addressing of the symbol coder ROM 21. The trellis decoder 15 
shown generally in FIGURE 1 is more particularly shown in FIGURE 2 as a 
Viterbi-type trellis decoder 40. 

Subsection 4.2.5 of ATSC Document A/53 Annex D I entitled 'Trellis coding" 
and prescribes a trellis coder for use in a DTV transmitter. The FIGURE 2 DTV receiver 
could be modified to use such a trellis coder for responding to the data that the Viterbi- 
type trellis decoder 40 generates for application to the data de-interleaver 16, the trellis 
coder generating a symbol stream for application to the interpolation filter 22 as its 
respective input signal. Rather than using a trellis coder to generate the symbol stream 
applied to the interpolation filter 22 as its respective input signal, however, FIGURE 2 
shows preliminary data-slicing results being extracted from the Viterbi-type trellis 
decoder 40 for application to the symbol coder ROM 21 as input addressing. The 
compensatory digital delay 24 is replaced by compensatory digital delay 41 providing 
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delay compensating for the combined latent delays of the adaptive filter 13, the 
decimation filter circuitry 14, the Viterbi-type trellis decoder 40 up to its preliminary 
data-slicing results output port, the symbol coder ROM 21, and the interpolation filter 
circuitry 22. 

5 FIGURE 4 diagrams a DTV receiver that differs from the FIGURE 1 DTV 

receiver in that an interpolation filter 42, which is an FIR digital filter with a system 
function that corresponds to the ideal channel system response, replaces the interpolation 
filter 22 with a system function that includes lowpass filtering for suppressing aliasing 
engendered by re-sampling to higher sample rate, but not otherwise affecting the data 

10 slicer 20 response. The DFT computer 23 used in the FIGURE 1 DTV receiver to 
compute the DFT of the interpolation filter 22 response is replaced in the FIGURE 4 
DTV receiver by a DFT computer 43 of similar construction, used to compute the DFT of 
the interpolation filter 42 response. The digital delay 24 of the FIGURE 1 DTV receiver 
is replaced in the FIGURE 4 DTV receiver by a compensatory digital delay 44 for 

1 5 delaying the adaptive filter 13 input signal so as to be in temporal alignment with the 
interpolation filter 42 response. That is, the digital delay 44 compensates for the 
combined latent delays of the adaptive filter 13, the decimation filter circuitry 14, the data 
slicer 20, the symbol coder ROM 21 and the interpolation filter 42. The DFT computer 
25 used in the FIGURE 1 DTV receiver to compute the DFT of the adaptive filter 13 

20 input signal as delayed by the digital delay 24 is replaced in the FIGURE 4 DTV receiver 
by a DFT computer 45 of similar construction, used to compute the DFT of adaptive filter 
13 input signal as delayed by the digital delay 44. The DFT combiner 26 of the FIGURE 
1 DTV receiver is replaced in the FIGURE 4 DTV receiver by a similar-DFT-terms 
divider bank 46 for dividing the DFT of the interpolation filter 42 response, as computed 

25 by the DFT computer 43, term-by-corresponding-term from the DFT of the digital delay 
44 response to the adaptive-filter 13 input signal, as computed by the DFT computer 45. 
These term-by-corresponding-term divisions by the similar-DFT-terms divider bank 46 
generate a DFT that defines in the frequency domain the ideal adaptive filter 13 kernel 
weighting coefficients in the time domain. 
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The similar-DFT-terms divider bank 46 supplies the DFT of the adaptive filter 13 
kernel weighting coefficients that it generates to the I-DFT computer 28. The I-DFT 
computer 28 generates the time-domain adaptive filter 13 kernel weighting coefficients 
corresponding to the DFT supplied by the similar-DFT-terms divider bank 46. In the 
5 FIGURE 4 DTV receiver, as in the FIGURE 1 DTV receiver, these adaptive filter 13 

kernel weighting coefficients in the time domain are supplied to respective ones of a bank 
29 of digital lowpass filters which suppress fluctuations in the values of the adaptive filter 
13 kernel weighting coefficients caused by noise. 

FIGURE 5 is a flow chart showing in more detail a method 50 comprising steps 

1 0 51 -57 by which the adaptive filter 13 kernel weighting coefficients can be calculated in 
the FIGURE 4 DTV receiver. Steps 51, 53, 56 and 57 of the method 50 respectively 
correspond to the steps 31, 33, 37 and 38 of the method 30. Initially in the method 50, 
the DFT of the delayed adaptive-filter 13 input signal supplied from the compensatory 
digital delay 44 is computed by the DFT computer 45 in the step 51 using a fast Fourier 

1 5 transform algorithm of the well-known Cooley-Tukey type, and the estimation of the 
actually transmitted signal as re-sampled and lowpass filtered with ideal channel system 
function by the interpolation filter 42 is computed by the DFT computer 43 in the step 52. 
Following the step 51, the DFT of the delayed adaptive-filter 13 input signal is 
continuously re-computed by the DFT computer 45 in a repeated updating step 53 using a 

20 Fourier-shift method such as described in U. S. patent No. 4,972,358. Following the step 
52 using a fast Fourier transform algorithm of the well-known Cooley-Tukey type for 
initial computation of the DFT of the interpolation filter 42 response, the DFT of the 
interpolation filter 42 response is continuously re-computed in a repeated updating step 
54 performed by the DFT computer 43 using the same Fourier-shift method that the DFT 

25 computer 45 uses in performing the repeated updating step 53. 

The similar-DFT-terms divider bank 46 performs the step 55 of the method 50. In 
the step 55, the DFT of the interpolation filter 42 response as computed by the steps 52 
and 54 is divided term-by-corresponding-term by the DFT of the delayed 
adaptive-filter 13 input signal as computed by the steps 51 and 53. These divisions in the 
30 step 55 generate a DFT that defines in the frequency domain the ideal adaptive filter 13 
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kernel weighting coefficients in the time domain. Divisors that are so small they would 
cause out-of-range quotients should be augmented in these division procedures. The 
divisions in step 55 can be performed logarithmically using tables stored in read-only 
memory, for example. 

5 The computer 28 computes the inverse-DFT of this DFT in step 56 to obtain 

weighting coefficients for the adaptive filter 13 kernel. During initialization or 
re-initialization, these weighting coefficients can be directly applied to the weighting 
coefficients registers) of the adaptive filter 13 to speed up initialization or 
re-initialization as much as possible. Generally, however, in the step 57 each of the 

10 weighting coefficients is averaged over some time for reducing the effects of noise before 
being supplied to the weighting coefficients registers) of the adaptive filter 13, this 
averaging being done by the bank 29 of digital lowpass filters. 

FIGURE 6 diagrams a DTV receiver that differs from the FIGURE 4 DTV 
receiver in that it does not use data slicing of the adaptive filtering response to generate, 

15 in accordance with the so-called "hard"-decision technique, the input addressing of the 
symbol coder ROM 21 used for estimating actually transmitted VSB DTV signal 
symbols. Instead, the FIGURE 6 DTV receiver employs Viterbi decoding of the adaptive 
filtering response response to generate, in accordance with the so-called "soff '-decision 
technique, the input addressing of the symbol coder ROM 21. The trellis decoder 15 

20 shown generally in FIGURE 4 is more particularly shown in FIGURE 6 as a 

Viterbi-type trellis decoder 40. Rather than using a trellis coder to generate the symbol 
stream applied to the interpolation filter 42 as its respective input signal, FIGURE 6 
shows preliminary data-slicing results being extracted from the Viterbi-type trellis 
decoder 40 for application to the symbol coder ROM 21 as input addressing. The 

25 compensatory digital delay 44 is replaced by compensatory digital delay 47 providing 
delay compensating for the combined latent delays of the adaptive filter 13, the 
decimation filter circuitry 14, the Viterbi-type trellis decoder 40 up to its preliminary 
data-slicing results output port, the symbol coder ROM 21 and the interpolation filter 
circuitry 42. 
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The symbol coder ROM 21 used in the DTV signal receivers of FIGURES 1 and 
4 includes circuitry for parsing data supplied from the data slicer 20 into input addressing 
for the read-only memory. Parsing can be done relying on hand-shake techniques with 
the data slicer 20 ? or by relying on a data count supplied from the receiver data 
5 synchronizing circuitry. The symbol coder ROM 21 used in the DTV signal receivers of 
FIGURES 3 and 6 includes circuitry for parsing data supplied from the Viterbi trellis 
decoder 40 into input addressing for the read-only memory. Parsing can be done relying 
on hand-shake techniques with the Viterbi trellis decoder 40, or by relying on a data 
count supplied from the receiver data synchronizing circuitry. 

10 The data synchronizing circuitry is not particularly material to the invention and 

so is omitted from the drawing to avoid clutter therein. The data synchronizing circuitry 
usually comprises counters with counts reset responsive to pseudorandom noise 
sequences in the data field synchronizing signal and possibly to data segment 
synchronizing signal These synchronizing signals are preferably drawn from the 

15 over-sampled adaptive filter 13 response, before its decimation. The synchronizing 
signal should be over-sampled to avoid the loss of important phase information. 

As noted in the BACKGROUND OF INVENTION, the ultimate limit on how 
rapidly the weighting coefficients in the kernel of an auto-regressive type of adaptive 
filtering can be re-initialized after loss of tracking is determined by the desideratum that 

20 re-initialization of the weighting coefficients should be deferred until the data on which 
their computation is based exit the kernel storage registers) of the adaptive filtering. A 
similar ultimate limit applies to the types of adaptive filtering shown in FIGURES 1, 3, 4 
and 6, since the DFT of the adaptive filtering response is necessary to the computation of 
the weighting coefficients in the kernel of the adaptive filtering. The adaptation of the 

25 types of adaptive filtering shown in FIGURES 1, 3, 4 and 6 is faster than auto-regressive 
types, however, owing to the computation of weighting coefficients that depend upon the 
adaptive filtering response taking much less time. This is also noted in the 
BACKGROUND OF INVENTION. 
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FIGURE 7 shows in detail one form of the demodulator and over-sampling 
analog-to-digital conversion circuitry 12 used in the DTV receiver of FIGURE 1, 3, 4 or 
6. A local oscillator 120 is subject to automatic frequency and phase control (AFPC) of 
its local oscillations, which are used in down-converting to baseband the amplified 
intermediate-frequency signal supplied from the DTV receiver front-end 11. The local 
oscillator 120 supplies its oscillations in 0° phasing to a mixer 121 for mixing with the 
amplified I-F signal the DTV receiver front-end 11 supplies to the mixer 121. The mixer 
121 output signal is applied as input signal to an analog lowpass filter 122. The lowpass 
filter 122 selectively responds to a baseband down-conversion component of the mixer 
121 output signal, but does not respond to the image up-conversion component of the 
mixer 121 output signal. An over-sampling analog-to-digital converter 123 digitizes the 
lowpass filter 122 response to generate digitized in-phase baseband signal supplied as 
input signal to the adaptive filter 13 for equalization and echo-suppression. 

The local oscillator 120 supplies its oscillations in 90° phasing to a mixer 124 for 
mixing with the amplified I-F signal the DTV receiver front-end 11 supplies to the mixer 

124. The mixer 124 output signal is applied as input signal to an analog lowpass filter 

125. The lowpass filter 125 selectively responds to a baseband down-conversion 
component of the mixer 124 output signal, but does not respond to the image 
up-conversion component of the mixer 124 output signal. An over-sampling 
analog-to-digital converter 126 digitizes the lowpass filter 125 response to generate 
digitized quadrature-phase baseband signal supplied as input signal to an adaptive filter 
58 for equalization and echo-suppression. The adaptive filter 13 and the adaptive filter 
58 have similar kernels, both receiving updated weighting coefficients from the bank 29 
of digital lowpass filters in response to the time-domain kernel weighting coefficients that 
the I-DFT computer 28 generates. The response of the adaptive filter 58 is equalized, 
echo-free, imaginary-only baseband DTV signal. A digital-to-analog converter 59 
converts this signal to analog form for analog lowpass filtering by an AFPC-loop lowpass 
filter 60. The lowpass filter 60 response is supplied to the local oscillator 120 as its 
automatic frequency and phase control (AFPC) signal. This closes the AFPC feedback 
loop. 
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While the FIGURE 7 arrangement uses duplicate adaptive filters 13 and 58, the 
tracking control of local oscillator phase reduces the need for over-sampling in the 
equalization and echo-suppression filtering. Alternatively, rather than using separate 
adaptive filters 13 and 58 for the baseband signals generated by in-phase and 
quadrature-phase demodulation, these baseband signals can be time-division-multiplexed 
for reflexing through a single dual-phase adaptive filter, and then demultiplexed for 
respective application to the decimation filter circuitry 14 and to the DAC 59. This 
alternative halves the number of digital multipliers required in the adaptive filtering for 
equalization and echo-suppression, but doubles the clocking rate required in the adaptive 
filtering, of course. There are also techniques known that can extract AFPC information 
from the adaptive filter 13 real-only response if at least 100% oversampling is used. 

The DTV receivers of FIGURES 1, 3, 4 and 6 each employ baseband equalization 
and, when the carrier phase tracking methods described above with reference to FIGURE 
7 are used, the adaptive filtering 13 is included in the carrier-phase-tracker AFPC loop. 
There is general knowledge in the equalizer art that it is preferable to exclude the 
equalization and echo-suppression filter from the AFPC loop used for phase tracking, 
particularly if the filter is designed for suppressing long-delayed echoes. There is also 
general knowledge that passband equalization is a good general approach to employ in 
order to be able to exclude the adaptive filter from the AFPC loop. 

FIGURE 8 shows portions of a novel digital television receiver that employs an 
adaptive filter for passband equalization and echo-suppression, which adaptive filter has 
its kernel weighting coefficients computed in accordance with the invention. The DTV 
receiver front-end 61 amplifies received DTV signal and converts it in frequency to a 
low-band intermediate-frequency band offset from zero-frequency by a 500 kHz or so, in 
which low-band phase-splitting of the passband signal by digital filtering is feasible. The 
down-conversion is performed so the carrier is at the higher-frequency end of this 
low-band I-F band. However, because the DFT extends to zero frequency, the offset 
should not be too much more than 500 kHz or so, or a fairly large number of the 
lower-frequency spectral bins in the DFT will be empty of any information useful in DFT 
computations of kernel weighting coefficients of the adaptive filter. An over-sampling 
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analog-to-digital converter 62 is connected for receiving the analog low-band 
intermediate-frequency signal supplied from the DTV receiver front-end 61 and digitizes 
that analog signal to generate digital input signal for an adaptive filter 63 for providing 
passband equalization and echo-suppression in its response to its digital input signal. 
This response is supplied to circuitry 64 for performing a complex digital synchrodyne to 
generate an in-phase baseband response I descriptive of real-only demodulated DTV 
signal and to generate an quadrature-phase baseband response Q descriptive of the 
imaginary-only Hilbert transform of the real-only demodulated DTV signal. This 
complex digital synchrodyne is supported by digital descriptions of local oscillations in 
0° and 90° phasings supplied to the circuitry 64 from a digital controlled oscillator (DCO) 
65. The quadrature-phase baseband response Q from the complex digital synchrodyne 
circuitry 64 is supplied to a digital lowpass filter 66, which supplies a response used by 
the DCO 65 for automatic frequency and phase control of the descriptions of local 
oscillations in 0° and 90° phasings the DCO 65 supplies to the circuitry 64. That is, the 
digital lowpass filter 66 is used as the principal filter in the AFPC loop for the complex 
digital synchrodyne. The in-phase baseband response I from the complex digital 
synchrodyne circuitry 64 is supplied to decimation filter circuitry 67 for resampling to 
baud rate. 

FIGURE 8 shows the decimation filter circuitry 67 connected for supplying its 
response to a trellis decoder 68, which is of the twelve-phase type conventionally 
employed in an ATSC VSB DTV receiver. The trellis coder 68 supplies symbol 
decoding results to further conventional portions of the ATSC VSB DTV receiver, which 
are of the sort previously described with regard to the FIGURE 1 ATSC VSB DTV 
receiver. The decimation filter circuitry 67 is further connected for supplying its 
response to a data sheer 70 that detects data slices by a "hard"-decision symbol decoding 
method. The data slicer 70 supplies its "hard"-decision symbol decoding results as input 
addressing to a read-only memory 71 that recodes these decoding results using ideal 
symbol code levels. The symbol coder ROM 71 output signal provides estimates of the 
symbols that were actually transmitted by the VSB DTV OTA broadcast transmitter. The 
estimates of actually transmitted symbols provided by the symbol coder ROM 71 are 
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supplied to interpolation filter circuity 72 for re-sampling back to the sampling rate used 
in the adaptive filter 63. The response of the interpolation filter circuitry 72 is supplied 
as modulating input signal to a balanced amplitude modulator 73, which generates a 
double-sideband amplitude-modulation (DSB AM) signal for application to an 
ideal-channel-response vestigial-sideband filter 74. The VSB filter 74 not only 
suppresses the higher-frequency sideband of the DSB AM signal it receives as input 
signal, but also suppresses the lowest-frequency portions of the remaining sideband to 
establish the Nyquist slope for reducing inter-symbol interference (ISI). The 
VSB filter 74 response estimates the VSB signal actually transmitted by the transmitter. 
DFT computation circuitry 75 computes the DFT of the estimated VSB signal actually 
transmitted by the transmitter. 

The balanced amplitude modulator 73 comprises a digital multiplier receiving, as 
its multiplicand and multiplier input signals, the symbol coder ROM 71 output signal and 
digitized carrier signal from the DCO 65. The phase of the digitized carrier signal is 
retarded from 0° to compensate for the combined latent delays of the adaptive filter 63, 
the digital synchrodyne circuitry 64, the decimation filter circuitry 67, the data slicer 70, 
the symbol coder ROM 71, the interpolation filter circuitry 72, the balanced amplitude 
modulator 73 and the VSB filter 74. The boundaries of the data slicing by the data slicer 
70 preferably take into account the pilot carrier accompanying the received symbols. 
That is, the boundaries of the data slicing by the data slicer 70 are asymmetrically 
disposed relative to the zero-energy axis for detecting normalized modulation levels of 
-5.75, -3.75, -1.75, +0.25, +2.25, +4.25. +6.25, + 8.25 nominal values. The symbol 
coder ROM 71 output signal augments the modulation level, as quantized to -5.75, -3.75, 
-1.75, +0.25, +2.25, +4.25. +6.25 and + 8.25 values, with an offset designed to offset the 
roll-off of pilot carrier component that will occur in the VSB filter 74. The asymmetry of 
the modulation levels of the symbol coder ROM 71 output signal with regard to zero 
causes a direct component in the system response of this signal applied to the balanced 
amplitude modulator 73, which unbalances the modulator to insert a fixed-amplitude pilot 
carrier into its digitized DSB AM output signal. The symbol coder ROM 71 in the 
FIGURE 8 DTV signal receiver includes circuitry for parsing data supplied from the data 
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slicer 70 into input addressing for the read-only memory. Parsing can be done relying on 
hand-shake techniques with the data slicer 70, or by relying on a data count supplied from 
the receiver data synchronizing circuitry. The synchrodyne circuitry 64 supplies its 
in-phase output signal I to the receiver date synchronizing circuitry as input signal 
5 thereto. 

The estimation of the signal that the VSB DTV OTA broadcast transmitter 
actually transmitted to the FIGURE 8 VSB DTV receiver, as supplied by the VSB filter 
74, is delayed respective to the actually received signal supplied from the over-sampling 
analog-to-digital conversion circuitry 62 to the adaptive filter 63 as input signal thereto. 

10 The input signal to the adaptive filter 63 is subjected to compensatory digital delay 76 
equal to the combined latent delays of the adaptive filter 63, the digital synchrodyne 
circuitry 64, the decimation filter circuitry 67, the data slicer 70, the symbol coder ROM 
71, the interpolation filter circuitry 72, the balanced amplitude modulator 73 and the VSB 
filter 74. DFT computation circuitry 77 computes the DFT of the delayed input signal to 

1 5 the adaptive filter 63 that the compensatory digital delay 76 provides. 

A similar-DFT-terms divider bank 78 is connected for dividing the DFT of the 
VSB filter 74 response, as computed by the DFT computer 75, 
term-by-corresponding-term by the DFT of the digital delay 76 response to the 
adaptive-filter 63 input signal, as computed by the DFT computer 77. These 
20 term-by-corresponding-term divisions by the similar-DFT-terms divider bank 78 generate 
a DFT that defines in the frequency domain the ideal adaptive filter 63 kernel weighting 
coefficients in the time domain. 

The similar-DFT-terms divider bank 78 supplies the DFT of the adaptive filter 63 
kernel weighting coefficients that it generates to an inverse discrete Fourier transform 
25 computer 79. The I-DFT computer 79 generates the time-domain adaptive filter 63 
kernel weighting coefficients corresponding to the DFT supplied by the 
similar-DFT-terms divider bank 78. These adaptive filter 63 kernel weighting 
coefficients in the time domain are supplied to respective ones of a bank 80 of digital 
lowpass filters that suppresses the effects of noise on the previous computation of the 
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kernel weighting coefficients. The bank 80 of digital lowpass filters supplies the 
lowpass-filtered adaptive filter 13 kernel weighting coefficients to coefficient registers in 
the adaptive filter 63 to be used as multiplier input signals by digital multipliers in the 
adaptive filter 63. The adaptive filter 63 is presumed to have the canonic transversal 
5 filter structure in which a tapped digital delay line supplies differently delayed tap 

responses to the adaptive filter 63 input signal for application to the digital multipliers as 
multiplicand input signals, and in which the products from the digital multipliers are 
summed to generate the adaptive filter 63 output signal. 

FIGURE 9 diagrams a DTV receiver that differs from the FIGURE 8 DTV 

1 0 receiver in that it does not use data slicing of the adaptive filtering response to generate, 
in accordance with the so-called "hard"-decision technique, the input addressing of the 
symbol coder ROM 71 used for estimating actually transmitted VSB DTV signal 
symbols. Instead, the FIGURE 9 DTV receiver employs Viterbi decoding of the adaptive 
filtering response response to generate, in accordance with the so-called "soft"-decision 

1 5 technique, the input addressing of the symbol coder ROM 71 . The trellis decoder 68 
shown generally in FIGURE 8 is more particularly shown in FIGURE 9 as a 
Viterbi-type trellis decoder 81. Rather than using a trellis coder to generate the symbol 
stream applied to the interpolation filter 72 as its respective input signal, FIGURE 9 
shows preliminary data-slicing results being extracted from the Viterbi-type trellis 

20 decoder 81 for application to the symbol coder ROM 71 as input addressing. The 
compensatory digital delay 76 is replaced by compensatory digital delay 82 providing 
delay compensating for the combined latent delays of the adaptive filter 63, the digital 
synchrodyne circuitry 64, the decimation filter circuitry 67, the Viterbi-type trellis 
decoder 81 up to its preliminary data-slicing results output port, the symbol coder ROM 

25 71 , the interpolation filter circuitry 72, the balanced amplitude modulator 73 and the VSB 
filter 74. The symbol coder ROM 71 used in the FIGURE 9 DTV signal receiver 
includes circuitry for parsing data supplied from the Viterbi trellis decoder 81 into input 
addressing for the read-only memory. Parsing can be done relying on hand-shake 
techniques with the Viterbi trellis decoder 81, or by relying on a data count supplied from 

30 the receiver data synchronizing circuitry. The synchrodyne circuitry 64 supplies its 
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in-phase output signal I to the receiver data synchronizing circuitry as input signal 
thereto. 

Variants of the FIGURE 8 and FIGURE 9 DTV signal receivers are contemplated 
in which the VSB filter 74 is not used, being replaced by direct connection or a wider 
band filter, and in which the ideal-channel-response VSB filtering is instead performed in 
the DFT regime. These variants include read-only memory for storing the DFT of the 
ideal-channel-response VSB filtering, and the similar-DFT-terms divider bank 78 is 
replaced by other similar-DFT-terms combining circuitry for combining the DFT of the 
ideal-channel-response VSB filtering with the DFTs computed by the DFT computers 75 
and 77. The combining procedures are in general the same as those used in the FIGURE 
1 and 3 DTV signal receivers. 

Thusfar in this specification, adaptation of the adaptive filtering used for 
equalization and echo-suppression has been based on real-only received signals. There 
are aspects of the invention that concern adaptive filtering used for equalization and 
echo-suppression of complex received signals, the adaptation of which adaptive filtering 
is based on discrete Fourier transforms of complex received signals. As an aid to 
understanding the DTV signal receiver apparatus described with reference to FIGURES 
10, 1 1, 12 and 13 it is useful to think of the transmitted VSB signal as comprising two 
component VSB signals modulating carriers at -45° and +45° phasing respective to the 0 C 
phasing of the actually transmitted VSB signal. Interestingly, this analysis results in two 
orthogonal channels with equal energies, rather than two orthogonal channels one with 
energy and the other without energy. The actually transmitted VSB signal is considered 
to be an in-phase carrier modulated by a real signal (Re) and to be accompanied by a 
zero- valued quadrature carrier modulated by an imaginary signal (Im). These two 
component VSB signals modulate the carriers at -45° and +45° phasing with modulating 
signals that have [(Re) 0 5 - j(Im)° 5 ] value and [(Re) 0 5 + jam) 0 5 ] value, respectively. 
Filtering the modulating signal that has [(Re) 0 5 + j(Im)° 5 ] value filters the real 
component (Re) component of complex received signal and the imaginary component 
(Im) of complex received signal similarly. Filtering the modulating signal that has 
[(Re) 0 5 - j(Im)° 5 ] value filters the real component (Re) and the imaginary component 
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(Im) differentially. If the two component filters implement equalization and 
echo-suppression correctly, the result of additively combining the responses of the two 
component filters generates a filter response for the real component (Re) of complex 
received signal. If the two component filters implement equalization and 

5 echo-suppression correctly, the result of subtractively combining the responses of the two 
component filters generates a zero-valued filter response for the imaginary component 
(Im) of complex received signal. If the two component filters implement equalization 
and echo-suppression correctly, the equalized signal as synchrodyned to baseband both at 
-45° phasing and at +45° phasing will in each case resemble the real-only signal (Re) as 

1 0 reduced in amplitude by square-root of two. This provides the basis for estimating the 
components of the actually transmitted signal at -45° phasing and at +45° phasing. This 
means that the two component filters can be adapted independently of each other. 

FIGURE 10 shows a DTV signal receiver that employs an adaptive filter for 
complex baseband equalization and echo-suppression, which adaptive filter has its kernel 

1 5 weighting coefficients computed in accordance with the invention. In the FIGURE 1 0 
DTV signal receiver the DTV receiver front-end 11 supplies amplified I-F signal to a 
mixer 83 for heterodyning with oscillations from a local oscillator 84. These oscillations 
are subject to automatic frequency and phase control (AFPC), their frequency being such 
that the mixer 83 output signal includes a low-band I-F signal. An analog lowpass filter 

20 85 is connected to selectively respond to this low-band I-F signal and not to its 

very-high-frequency image. (In a variant of the FIGURE 10 DTV signal receiver an 
analog bandpass filter that selectively responds to the low-band I-F signal replaces the 
lowpass filter 85.) The response of the filter 85 is applied as input signal to 
analog-to-digital conversion circuitry 86 for digitization. 

25 The ADC circuitry 86 is connected for supplying digitized low-band I-F signal to 

digital synchrodyne-to-baseband circuitry 87 arranged to synchronously detect the 
digitized low-band I-F signal at -45° phasing and at +45° phasing. The complete 
adaptive filtering used for equalization and echo suppression comprises: a component 
adaptive filter 88 connected to receive as its input signal the baseband signal the 
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synchrodyne circuitry 87 supplies as the result of synchronously detecting at -45 
phasing with respect to VSB DTV signal carrier, a component adaptive filter 89 
connected to receive as its input signal the baseband signal the synchrodyne circuitry 87 
supplies as the result of synchronously detecting at +45° phasing with respect to VSB 
5 DTV signal carrier, a digital adder 90 connected to sum the responses of the component 
adaptive filters 88 and 89 to generate a substantially-echo-free real equalized baseband 
DTV signal, and a digital subtracter 91 connected to differentially combine the responses 
of the component adaptive filters 88 and 89 to generate a substantially-echo-free 
imaginary equalized baseband DTV signal. 

10 The substantially-echo-free imaginary equalized baseband DTV signal the digital 

subtracter 91 supplies as its digital difference output signal should be zero-valued if the 
oscillations supplied by the local oscillator 84 are of correct frequency and phase, so that 
the digital synchrodyne-to-baseband circuitry 87 synchrodynes the digitized 
low-band I-F signal to baseband at -45° phasing and at +45° phasing. Otherwise, the 

1 5 substantially-echo-free imaginary equalized baseband DTV signal that the digital 

subtracter 91 supplies as difference output signal is a digital error signal for the AFPC 
loop that controls the frequency and phase of the oscillations the local oscillator 84 
supplies. The digital difference output signal from the digital subtracter 91 is supplied to 
a digital-to-analog converter 92 for conversion to an analog error signal applied to an 

20 AFPC-loop analog lowpass filter 93 as input signal. The response of the AFPC-loop 
lowpass filter 93 is applied to the local oscillator 84 as AFPC signal, thereby closing the 
degenerative-feedback AFPC loop that adjusts the frequency and phase of the oscillations 
the local oscillator 84 supplies. This adjustment is made so that the digital 
synchrodyne-to-baseband circuitry 87 synchrodynes the digitized low-band I-F signal to 

25 baseband at -45° phasing and at +45° phasing. Customary variations of the AFPC loop 
design, such as performing part of the lowpass AFPC-loop filtering in the digital domain, 
are possible, of course. 

The substantially-echo-free real equalized baseband DTV signal the digital adder 
90 supplies as its digital sum output signal is supplied to data synchronizing circuitry 94 
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and to the trellis decoder 15. The elements 15 - 19 in the FIGURE 10 DTV signal 
receiver correspond both in form and in function to the elements 15 - 19 in the FIGURE 
1 DTV signal receiver. 

A computer 95 computes the kernel weighting coefficients for the component 
adaptive filter 88, proceeding from the input signal supplied to the component adaptive 
filter 88 and from the response of the component adaptive filter 88. A computer 96 
computes the kernel weighting coefficients for the component adaptive filter 89, 
proceeding from the input signal supplied to the component adaptive filter 89 and from 
the response of the component adaptive filter 89. The computers 95 and 96 are alike in 
their construction, which is generally similar to that of computers used in the FIGURE 1, 
3, 4, and 6 DTV signal receivers for computing the kernel weighting coefficients for 
adaptive filtering of real-only signals. The difference in construction arises from the 
sampling rate in the component adaptive filters 88 and 89 presumably being baud rate, 
with complex sampling providing the "over-sampling" necessary to characterize the 
received signal as phase-modulated as well as amplitude-modulated by multipath 
distortion. Accordingly, there is no decimation filter for re-sampling adaptive filter 
response before the data slicing or Viterbi decoding, and there is no interpolation filter 
for re-sampling the baud-rate symbol decoder ROM response back to a higher rate of 
sampling used in the adaptive filter. Of course, complex adaptive filtering designs that 
use higher-than-baud-rate sampling can be constructed using the principles of the 
invention, in which case the computers 95 and 96 will incorporate decimation filters 
before symbol decoding and interpolation filters after re-coding symbols. 

FIGURE 1 1 shows a further DTV signal receiver that employs an adaptive filter 
for complex baseband equalization and echo-suppression, which adaptive filter has its 
kernel weighting coefficients computed in accordance with the invention. The FIGURE 
1 1 DTV signal receiver differs from the FIGURE 10 DTV signal receiver in that the 
DTV receiver front-end 61 replaces the DTV receiver front-end 11, the mixer 83, the 
AFPC'd local oscillator 84, and the analog lowpass filter 85. The DTV receiver 
front-end 61 converts the DTV signal to a low-band intermediate-frequency band offset 
from zero-frequency by a 500 kHz or so. The FIGURE 1 1 DTV signal receiver also 
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differs from the FIGURE 10 DTV signal receiver in that the FIGURE 10 digital 
synchrodyne-to-baseband circuitry 87 that synchrodynes to baseband using orthogonal 
digital carriers with fixed-frequency-and-phase system characteristics is replaced by- 
digital synchrodyne-to-baseband circuitry 97 that synchrodynes to baseband using 

5 orthogonal digital carriers with controlled-frequency-and-phase system characteristics as 
supplied from a digital controlled oscillator 98. The FIGURE 1 1 DTV signal receiver 
dispenses with the DAC 92 and the AFPC-loop analog lowpass filter 93 used in the 
FIGURE 10 DTV signal receiver to complete the AFPC loop for the controlled local 
oscillator 84. Instead, the AFPC loop of the FIGURE 1 1 DTV signal receiver uses an 

10 AFPC-loop digital lowpass filter 99 to complete an AFPC loop for the DCO 98. More 
particularly, the digital difference output signal from the digital subtracter 91 is supplied 
to the digital lowpass filter 99 as input signal, and the response of the lowpass filter 99 is 
applied to the DCO 98 as AFPC signal. This closes the degenerative-feedback AFPC 
loop that adjusts the frequency and phase of the oscillations the DCO 98 supplies to the 

1 5 digital synchrodyne-to-baseband circuitry 97 so that circuitry 97 synchrodynes the 

digitized low-band I-F signal from the ADC 86 to baseband at -45° phasing and at +45° 
phasing. 

The DTV receivers of FIGURES 10 and 1 1 each employ baseband equalization 
and the complex adaptive filtering provided by the elements 88 - 91 is included in the 
20 carrier-phase-tracker AFPC loop. As remarked earlier in the specification, it is preferable 
to exclude the equalization and echo-suppression filter from the AFPC loop used for 
phase tracking, and passband equalization is a good general approach to employ in order 
to be able to exclude the adaptive filter from the AFPC loop. 

FIGURE 12 shows portions of a novel digital television receiver that employs a 
25 complex adaptive filter for passband equalization and echo-suppression, which adaptive 
filter has its kernel weighting coefficients computed in accordance with the invention. In 
the FIGURE 12 DTV signal receiver, just as in the FIGURE 10 DTV signal receiver, the 
DTV receiver front-end 11 supplies amplified I-F signal to a mixer 83 for heterodyning 
with oscillations from the controlled local oscillator 84. The analog lowpass filter 85 (or 



38 



Atty. Dkt. ALRL03A 



PATENT Application of J. D. McDonald et al. 



an analog bandpass filter) is connected for selecting the low-band I-F signal component 
of the mixer 83 output signal to the ADC circuitry 86 for digitization. 

In the FIGURE 12 DTV signal receiver the digitized low-band I-F signal is not 
directly supplied to the digital synchrodyne-to-baseband circuitry for demodulation 

5 before equalization and echo-suppression, in which digital synchrodyne-to-baseband 
circuitry the digitized low-band I-F signal is phase-split in preparation for the actual 
synchrodyning to baseband procedures. In the FIGURE 12 DTV signal receiver the ADC 
circuitry 86 supplies the digitized low-band I-F signal as input signal to a phase-splitter 
100. The phase-splitter 100 supplies orthogonal responses to the digitized low-band I-F 

10 signal as input signals to an adaptive filter 101 and to an adaptive filter 102. The 
responses of the adaptive filters 101 and 102 are supplied to digital 
synchrodyne-to-baseband circuitry 103 that synchrodynes these responses to baseband in 
orthogonal phasing to generate demodulation results from demodulating nominally at 
-45° phasing and at +45° phasing, respectively. 

1 5 A digital subtracter 104 differentially combines these baseband responses to 

generate an imaginary response for the complete adaptive filter for equalization and 
echo-suppression. This imaginary response is used as a digital error signal indicative of 
the corrections that need to be made in the frequency and phase of the local oscillations 
that the controlled local oscillator 84 supplies. Digital-to-analog conversion circuitry 105 

20 converts this digital error signal to an analog error signal applied as input signal to the 
AFPC-loop analog lowpass filter 93 as input signal. The response of the AFPC-loop 
lowpass filter 93 is applied to the local oscillator 84 as AFPC signal, thereby closing the 
degenerative-feedback AFPC loop that adjusts the frequency and phase of the oscillations 
the local oscillator 84 supplies. This adjustment is made so that the digital 

25 synchrodyne-to-baseband circuitry 103 synchrodynes the digitized low-band I-F signal to 
baseband by demodulating at -45° phasing and at +45° phasing. As in the FIGURE 10 
DTV signal receiver, variations of the AFPC loop design are possible, per custom. 

The digital synchrodyne-to-baseband circuitry 103 supplies its baseband 
responses obtained by demodulating at -45° phasing and at +45° phasing as input signals 
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for decimation filters 106 and 107, respectively, which re-sample these signals to baud 
rate. A digital adder 108 sums these re-sampled baseband responses to generate a real 
response for the complete adaptive filter for equalization and echo-suppression. The 
substantially-echo-free real equalized baseband DTV signal the digital adder 108 supplies 
as its digital sum output signal is supplied to data synchronizing circuitry 94 and to the 
trellis decoder 15. The elements 15 - 19 in the FIGURE 12 DTV signal receiver 
correspond both in form and in function to the elements 15 - 19 in the FIGURE 1 and 
FIGURE 10 DTV signal receivers. 

The kernel weighting coefficients for the component adaptive filter 101 are 
computed by a computer 109, responsive to component adaptive filter 101 input signal 
and to estimation of the actually transmitted signal giving rise to the component adaptive 
filter 101 input signal, which estimation is based on data-slicer response or Viterbi 
decoder response to the decimation filter 106 response. The kernel weighting 
coefficients for the component adaptive filter 102 are computed by a computer 110, 
responsive to component adaptive filter 102 input signal and to estimation of the actually 
transmitted signal giving rise to the component adaptive filter 102 input signal, which 
estimation is based on data-slicer response or Viterbi decoder response to the decimation 
filter 107 response. More particularly, the computation of the kernel weighting 
coefficients for each of the component adaptive filters 101 and 102 is generally similar to 
the computation of the kernel weighting coefficients of the adaptive filter 63 in the 
FIGURE 8 and FIGURE 9 DTV signal receivers. 

FIGURE 13 shows a further DTV signal receiver that employs an adaptive filter 
for complex passband equalization and echo-suppression, which adaptive filter has its 
kernel weighting coefficients computed in accordance with the invention. The FIGURE 
13 DTV signal receiver differs from the FIGURE 12 DTV signal receiver in that the 
DTV receiver front-end 61 replaces the DTV receiver front-end 11, the mixer 83, the 
AFPC'd local oscillator 84, and the analog lowpass filter 85. The DTV receiver 
front-end 61 converts the DTV signal to a low-band intermediate-frequency band offset 
from zero-frequency by a 500 kHz or so. The FIGURE 13 DTV signal receiver also 
differs from the FIGURE 12 DTV signal receiver in that the FIGURE 12 digital 
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synchrodyne-to-baseband circuitry 103 that synchrodynes to baseband using orthogonal 
digital carriers with fixed-frequency-and-phase system characteristics is replaced by 
digital synchrodyne-to-baseband circuitry 113 that synchrodynes to baseband using 
orthogonal digital carriers with controlled-frequency-and-phase system characteristics as 
supplied from a digital controlled oscillator 114. The FIGURE 13 DTV signal receiver 
dispenses with the DAC 106 and the AFPC-loop analog lowpass filter 93 used in the 
FIGURE 12 DTV signal receiver to complete the AFPC loop for the controlled local 
oscillator 84. Instead, the AFPC loop of the FIGURE 13 DTV signal receiver uses an 
AFPC-loop digital lowpass filter 115 to complete an AFPC loop for the DCO 114. More 
particularly, the digital difference output signal from the digital subtracter 104 is supplied 
to the digital lowpass filter 115 as input signal, and the response of the lowpass filter 115 
is applied to the DCO 114 as AFPC signal. This closes the degenerative-feedback AFPC 
loop that adjusts the frequency and phase of the oscillations the DCO 114 supplies to the 
digital synchrodyne-to-baseband circuitry 113 so that circuitry 113 synchrodynes the 
response of the adaptive filter 101 to baseband at -45° phasing and synchrodynes the 
response of the adaptive filter 102 to baseband at +45° phasing. 

FIGURE 14 shows a modification that is made in the adaptive filtering in the 
FIGURE 1, 3, 4 and 6 DTV signal receivers in other embodiments of the invention. In 
the modified FIGURE 1 and 3 DTV signal receivers, the cascade arrangement of 
inverse-DFT computer 28 responsive to the DFT generated by the similar-DFT-terms 
combiner 26 and the bank 29 of respective digital lowpass filters for the weighting 
coefficients generated by the computer 28 are replaced by a bank 116 of respective digital 
lowpass filters for the terms of the DFT generated by the similar-DFT-terms combiner 26 
and by an inverse-DFT computer 117 for generating the kernel weighting coefficients of 
the adaptive filter 13 from the responses of the bank 116 of respective digital lowpass 
filters. In the modified FIGURE 4 and 6 DTV signal receivers, the cascade arrangement 
of inverse-DFT computer 28 responsive to the DFT generated by the similar-DFT-terms 
divider bank 46 and the bank 29 of respective digital lowpass filters for the weighting 
coefficients generated by the computer 28 are replaced by a bank 118 of respective digital 
lowpass filters for the terms of the DFT generated by the similar-DFT-terms divider bank 
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46 and by an inverse-DFT computer 119 for generating the kernel weighting coefficients 
of the adaptive filter 13 from the responses of the bank 118 of respective digital lowpass 
filters. 

FIGURE 1 5 shows a modification that is made in the adaptive filtering in the 
FIGURE 8 and 9 DTV signal receivers in other embodiments of the invention. In the 
modified FIGURE 8 and 9 DTV signal receivers, the cascade arrangement of 
inverse-DFT computer 79 responsive to the DFT generated by the similar-DFT-terms 
divider bank 78 and the bank 80 of respective digital lowpass filters for the weighting 
coefficients generated by the computer 79 are replaced by a bank 118 of respective digital 
lowpass filters for the terms of the DFT generated by the similar-DFT-terms divider bank 
78 and by an inverse-DFT computer 119 for generating the kernel weighting coefficients 
of the adaptive filter 63 from the responses of the bank 118 of respective digital lowpass 
filters. 

In the foregoing description of the invention it was presumed that the adaptive 
digital filtering used for channel equalization and echo suppression was of FIR type. A 
problem that is encountered when using FIR filtering for echo suppression is that the 
signals with echoes from times previous or subsequent to the current time are weighted 
and combined with the current signal in order to suppress the echoes attending it. The 
desideratum for echo suppression is that echo-free signals from times previous or 
subsequent to the current time are weighted and combined with the current signal in order 
to suppress the echoes attending it. The unwanted weighted echo components in the 
combining procedure create repeats of the original echo. These repeats are an increasing 
problem as the echo components increase in amplitude vis-a-vis the principal signal. A 
procedure that can be done to reduce the errors in the kernel of the FIR filtering is as 
follows. The time-domain response of the actual reception channel to unit impulse is 
modified by nulling the principal signal term. This modified time-domain response is 
convolved with the zero-order set of originally determined weighting coefficients to 
develop a first-order set of correction weighting coefficients to be combined with the 
zero-order set of original weighting coefficients for canceling terms of that zero-order set. 
The modified time-domain response is convolved with the first-order set of correction 
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weighting coefficients to develop a second-order set of correction weighting coefficients 
to be combined with the zero-order set of original weighting coefficients and the 
first-order set of correction weighting coefficients for canceling terms of that first-order 
set. This successive correction procedure continues until additional correction weighting 
coefficients fall below the quantizing noise of the weighting coefficients. These 
procedures are best implemented performing the convolution procedures in the DFT 
domain. 

Some of the problem that is encountered when using FIR filtering for echo 
suppression is avoided in adaptive digital filtering of different structure that is more 
commonly used for channel equalization and echo suppression in DTV receivers. This 
alternative structure cascades a first FIR adaptive digital filter with an IIR adaptive digital 
filter, which is composed of a digital adder and a second FIR adaptive digital filter. The 
adder is connected to receive the IIR filter input signal as its first addend input signal and 
to supply the IIR filter output signal as its sum output signal. The second FIR digital 
filter is connected for feeding back the IIR filter output signal for combination with the 
IIR filter input signal. More particularly, the second FIR digital filter is connected to 
receive the IIR filter output signal as its input signal and to supply its response to the 
adder as a second addend input signal thereof. The IIR filter is used to suppress 
post-echoes, especially those with substantial delay. The first FIR filter is used to 
suppress pre-echoes and post-echoes with too short a delay to be conveniently suppressed 
by the IIR filter. The advantage in using the IIR filter to suppress post-echoes is that 
spectral repeats are automatically suppressed, so many fewer non-zero weighting 
coefficients are required for suppressing post-echoes than if completely FIR filtering 
were used instead. 

The weighting coefficients for the first and second FIR adaptive digital filters in 
this type of adaptive filtering can be computed by modifying the methods described 
above. The weighting coefficients associated with the suppression of post-echoes 
determined in the final step of the FIGURE 2 method or the FIGURE 5 method are 
supplied to the coefficient register of the second FIR filter. The zero-order set of 
originally determined weighting coefficients in the kernel of the first FIR filter comprises 
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the weighting coefficients associated with the suppression of pre-echoes and the 
weighting coefficient associated with the principal signal. 

The procedure that can then be done to reduce the errors in the kernel of the first 
FIR filter is as follows. The time-domain response of the actual reception channel to unit 

5 impulse is modified by nulling the principal signal term and each post-echo term. The 
resulting time-domain pre-echo-only response is convolved with the zero-order set of 
originally determined weighting coefficients for the first FIR filter to develop a first-order 
set of correction weighting coefficients for the first FIR filter to be combined with its 
zero-order set of original weighting coefficients for canceling terms of that zero-order set. 

1 0 The time-domain pre-echo-only response is convolved with the first-order set of 

correction weighting coefficients for the first FIR filter to develop a second-order set of 
correction weighting coefficients for the first FIR filter to be combined with its zero- 
order set of original weighting coefficients and its first-order set of correction weighting 
coefficients for canceling terms of that first-order set. This successive correction 

1 5 procedure continues until additional correction weighting coefficients fall below the 
quantizing noise of the weighting coefficients. These procedures are best implemented 
performing the convolution procedures in the DFT domain. 

In still other embodiments of the invention, the DFT computers shown explicitly 
in FIGURES 1, 3, 4, 6, 8 and 9 (and implicitly in FIGURES 10-13) are constructed as 

20 parallel-bank analyzer filters; and the inverse DFT computers shown explicitly in 
FIGURES 1, 3, 4, 6, 7, 8, 9, 14 and 15 (and implicitly in FIGURES 10-13) are 
constructed as parallel-bank synthesizer filters. Methods for constructing the analyzer 
and synthesizer filters are described at length in the textbook Multirate Digital Signal 
Processing, published in 1983 by Prentice-Hall, Inc. of Englewood Cliffs, NJ, and 

25 authored by Ronald E. Crochiere and Lawrence R. Rabiner. See particularly Chapter 7 
entitled "Multirate Techniques in Filter Banks and Spectrum Analyzers and Synthesizers' 
running from page 289 to page 404. 

Each of the FIGURES 16 - 19 show an adaptive filter structure 130 using such 
analyzer and synthesizer filters. The transversal adaptive filter 13 is replaced by an 
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analyzer filter 131 that converts the adaptive filter input signal to an analyzer filter 
response that is the DFT of the adaptive filter input signal, similar-DFT-terms multiplier 
bank 132 that multiplies the DFT of the adaptive filter input signal 
term-by-corresponding-term with the DFT of the desired system characteristic of the 
adaptive filter to generate the DFT of the adaptive filter response, and a synthesizer filter 
133 receptive of the resulting DFT of the adaptive filter response to generate the adaptive 
filter response in the time domain. 

FIGURE 16 shows how the adaptive filter structure 130 replaces the transversal 
adaptive filter 13 in a modification of the FIGURE 1 DTV signal receiver. The DFT of 
the desired system characteristic of the adaptive filter structure 130 is generated by 
methods similar to those used to generate the DFT of the kernel of the transversal 
adaptive filter structure in FIGURE 1 . The DFT computer 23 of the FIGURE 1 DTV 
signal receiver is replaced in the FIGURE 16 DTV signal receiver by an analyzer filter 
023 for generating the DFT of the interpolation filter 22 response to the estimates of 
transmitted symbols supplied by the symbol coder ROM 21. In the FIGURE 16 DTV 
signal receiver, the analyzer filter 131 and a parallel digital-delay bank 134 providing 
equal digital delays for each DFT term from the analyzer filter 131 replaces the digital 
delay 24 and the DFT computer 25 of the FIGURE DTV 1 signal receiver. The parallel 
digital-delay bank 134 compensates for the combined latencies of the synthesizer filter 
133, the decimation filter 14, the data slicer 20, the symbol coder ROM 21 and the 
interpolation filter 22 (as well as any difference between the latencies of the analyzer 
filters 131 and 023, should such difference exist). The inverse-DFT computer 28 and the 
bank 29 of digital lowpass filters used in the FIGURE 1 DTV signal receiver are 
dispensed with in the FIGURE 16 DTV signal receiver. Instead, the DFT of the system 
response for the adaptive filter structure 130 is generated by the bank 116 of digital 
lowpass filters responding to the terms of the DFT generated by the similar-DFT-terms 
combiner 26. This DFT of the adaptive filter structure 130 system response is supplied to 
the similar-DFT-terms multiplier bank 132 to be multiplied term-by-corresponding term 
by the DFT of the adaptive filter structure 130 input signal supplied by the analyzer filter 
131. The synthesizer filter 133 responds to the product output signal from the 
similar-DFT-terms multiplier bank 132 to generate the substantially-echo-free, equalized 
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output signal of the adaptive filter structure 130, which is supplied to the trellis decoder 
15 for symbol decoding and to the data sheer 20. The data sheer 20 supplies addressing 
to the symbol coder ROM 21, which ROM 21 generates estimates of the transmitted 
symbol code at baud rate. The interpolation filter 22 responds to these baud-rate 
estimates of the transmitted symbol code, re-sampling these estimates to the sample rate 
of the adaptive filter structure 130 input signal and supplying them to the analyzer filter 
023 as input signal thereto. 

FIGURE 17 shows how the adaptive filter structure 130 replaces the transversal 
adaptive filter 13 in a modification of the FIGURE 3 DTV signal receiver that is a still 
further embodiment of the invention. The DFT of the desired system characteristic of the 
adaptive filter structure 130 is generated by methods similar to those used to generate the 
DFT of the kernel of the transversal adaptive filter structure in FIGURE 3. The DFT 
computer 23 of the FIGURE 3 DTV signal receiver is replaced in the FIGURE 16 DTV 
signal receiver by an analyzer filter 023 for generating the DFT of the interpolation filter 
22 response to the estimates of transmitted symbols supplied by the symbol coder ROM 
21. In the FIGURE 17 DTV signal receiver, the analyzer filter 131 and a parallel 
digital-delay bank 135 providing equal digital delays for each DFT term from the 
analyzer filter 131 replaces the digital delay 41 and the DFT computer 25 of the FIGURE 
DTV 3 signal receiver. The parallel digital-delay bank 135 compensates for the 
combined latencies of the synthesizer filter 133, the decimation filter 14, the Viterbi 
trellis decoder 40, the symbol coder ROM 21 and the interpolation filter 22 (as well as 
any difference between the latencies of the analyzer filters 131 and 023, should such 
difference exist). The inverse-DFT computer 28 and the bank 29 of digital lowpass 
filters used in the FIGURE 3 DTV signal receiver are dispensed with in the FIGURE 17 
DTV signal receiver. Instead, the DFT of the system response for the adaptive filter 
structure 130 is generated by the bank 116 of digital lowpass filters responding to the 
terms of the DFT generated by the similar-DFT-terms combiner 26. This DFT of the 
adaptive filter structure 130 system response is supplied to the similar-DFT-terms 
multiplier bank 132 to be multiplied term-by-corresponding term by the DFT of the 
adaptive filter structure 130 input signal as supplied by the analyzer filter 131. The 
synthesizer filter 133 responds to the product output signal from the similar-DFT-terms 
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multiplier bank 132 to generate the substantially-echo-free, equalized output signal of the 
adaptive filter structure 130, which is supplied to the Viterbi trellis decoder 40 for symbol 
decoding. 

FIGURE 18 shows how the adaptive filter structure 130 replaces the transversal 
adaptive filter 13 in a modification of the FIGURE 4 DTV signal receiver. The DFT of 
the desired system characteristic of the adaptive filter structure 130 is generated by 
methods similar to those used to generate the DFT of the kernel of the transversal 
adaptive filter structure in FIGURE 4. The DFT computer 43 of the FIGURE 4 DTV 
signal receiver is replaced in the FIGURE 18 DTV signal receiver by an analyzer filter 
043 for generating the DFT of the interpolation filter 42 ideal-channel-system-function 
response to the estimates of transmitted symbols supplied by the symbol coder ROM 21. 
In the FIGURE 18 DTV signal receiver, the analyzer filter 131 and a parallel 
digital-delay bank 136 providing equal digital delays for each DFT term from the 
analyzer filter 131 replaces the digital delay 44 and the DFT computer 45 of the FIGURE 
DTV 4 signal receiver. The parallel digital-delay bank 136 compensates for the 
combined latencies of the synthesizer filter 133, the decimation filter 14, the data sheer 
20, the symbol coder ROM 21 and the interpolation filter 42 (as well as any difference 
between the latencies of the analyzer filters 131 and 023, should such difference exist). 
The inverse-DFT computer 28 and the bank 29 of digital lowpass filters used in the 
FIGURE 4 DTV signal receiver are dispensed with in the FIGURE 18 DTV signal 
receiver. Instead, the DFT of the system response for the adaptive filter structure 130 is 
generated by the bank 116 of digital lowpass filters responding to the terms of the DFT 
generated by the similar-DFT-terms divider bank 46. This DFT of the system response 
for the adaptive filter structure 130 is supplied to the similar-DFT-terms multiplier bank 
132 to be multiplied term-by-corresponding term by the DFT of adaptive filter structure 
130 input signal supplied by the analyzer filter 131. The synthesizer filter 133 responds 
to the product output signal from the similar-DFT-terms multiplier bank 132 to generate 
the substantially-echo-free, equalized output signal of the adaptive filter structure 130, 
which is supplied to the trellis decoder 15 for symbol decoding and to the data slicer 20. 
The data slicer 20 supplies addressing to the symbol coder ROM 21, which ROM 21 
generates estimates of the transmitted symbol code at baud rate. The interpolation filter 
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42 responds to these baud-rate estimates of the transmitted symbol code, re-sampling 
these estimates to the sample rate of the adaptive filter structure 130 input signal and 
supplying them to the analyzer filter 043 as input signal thereto. 

FIGURE 19 shows how the adaptive filter structure 130 replaces the transversal 
adaptive filter 13 in a modification of the FIGURE 6 DTV signal receiver that is a still 
further embodiment of the invention. The DFT of the desired system characteristic of the 
adaptive filter structure 130 is generated by methods similar to those used to generate the 
DFT of the kernel of the transversal adaptive filter structure in FIGURE 6. The DFT 
computer 43 of the FIGURE 6 DTV signal receiver is replaced in the FIGURE 16 DTV 
signal receiver by an analyzer filter 043 for generating the DFT of the interpolation filter 
42 ideal-channel-system-function response to the estimates of transmitted symbols 
supplied by the symbol coder ROM 21. In the FIGURE 19 DTV signal receiver, the 
analyzer filter 131 and a parallel digital-delay bank 137 providing equal digital delays for 
each DFT term from the analyzer filter 131 replaces the digital delay 47 and the DFT 
computer 45 of the FIGURE DTV 6 signal receiver. The parallel digital-delay bank 137 
compensates for the combined latencies of the synthesizer filter 133, the decimation filter 
14, the Viterbi trellis decoder 40, the symbol coder ROM 21 and the interpolation filter 
42 (as well as any difference between the latencies of the analyzer filters 131 and 023, 
should such difference exist). The inverse-DFT computer 28 and the bank 29 of digital 
lowpass filters used in the FIGURE 6 DTV signal receiver are dispensed with in the 
FIGURE 19 DTV signal receiver. Instead, the DFT of the system response for the 
adaptive filter structure 130 is generated by the bank 116 of digital lowpass filters 
responding to the terms of the DFT generated by the similar-DFT-terms divider bank 46. 
This DFT of the adaptive filter structure 130 system response is supplied to the 
similar-DFT-terms multiplier bank 132 to be multiplied term-by-corresponding term by 
the DFT of the adaptive filter structure 130 input signal supplied by the analyzer filter 
131. The synthesizer filter 133 responds to the product output signal from the 
similar-DFT-terms multiplier bank 132 to generate the substantially-echo-free, equalized 
output signal of the adaptive filter structure 130, which is supplied to the Viterbi trellis 
decoder 40 for symbol decoding. 
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As noted above, when adaptive FIR filtering is used for echo suppression and 
channel equalization, a problem arises because signals with echoes from times previous 
or subsequent to the current time are weighted and combined with the current signal in 
order to suppress the echoes attending it, rather than echo-free signals being so weighted 
and combined. As noted above, further convolution procedures carried out in the DFT 
domain can be used to generate corrections in the weighting coefficients computed for 
the kernel of the adaptive FIR filtering, so performance is improved especially when 
multipath distortion is severe. An analogous problem arises in regard to the weighting 
coefficients supplied to the bank 116 of digital lowpass filters in FIGURES 16, 17, 18 
and 19. The solution to this analogous problem is provided by further convolution 
procedures carried out in the DFT domain, which further procedures are analogous to 
those used for correcting the weighting coefficients computed for the kernel of the 
adaptive FIR filtering. 

FIGURE 20 shows a modification of the FIGURE 8 DTV signal receiver to 
replace the transversal adaptive filter 63 by an alternative adaptive filter structure 630. 
The adaptive filter structure 630 comprises an analyzer filter 631 that converts the 
adaptive filter input signal to an analyzer filter response that is the DFT of the adaptive 
filter input signal, similar-DFT-terms multiplier bank 632 that multiplies the DFT of the 
adaptive filter input signal term-by-corresponding-term with the DFT of the desired 
system characteristic of the adaptive filter to generate the DFT of the adaptive filter 
response, and a synthesizer filter 633 receptive of the resulting DFT of the adaptive filter 
response to generate the adaptive filter response in the time domain. The DFT of the 
desired system characteristic of the adaptive filter structure 630 is generated by methods 
similar to those used to generate the DFT of the kernel of the transversal adaptive filter 
structure in FIGURE 8. The DFT computer 75 of the FIGURE 8 DTV signal receiver is 
replaced in the FIGURE 20 DTV signal receiver by an analyzer filter 075 for generating 
the DFT of the ideal-channel-response VSB filter 74 output signal. In the FIGURE 20 
DTV signal receiver, the analyzer filter 631 and a parallel digital-delay bank 634 
providing equal digital delays for each DFT term from the analyzer filter 631 replaces the 
digital delay 76 and the DFT computer 77 of the FIGURE DTV 8 signal receiver. The 
parallel digital-delay bank 634 compensates for the combined latencies of the synthesizer 
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filter 633, the digital synchrodyne to baseband circuitry 64, the decimation filter 67, the 
data slicer 20, the symbol coder ROM 21,the interpolation filter 72, the balanced 
amplitude modulator 73, and the ideal-channel-response VSB filter 74 (as well as any 
difference between the latencies of the analyzer filters 631 and 075, should such 
difference exist). The inverse-DFT computer 79 and the bank 80 of digital iowpass 
filters used in the FIGURE 8 DTV signal receiver are dispensed with in the FIGURE 20 
DTV signal receiver. Instead, the DFT of the system response for the adaptive filter 
structure 630 is generated by the bank 118 of digital Iowpass filters responding to the 
terms of the DFT generated by the similar-DFT-terms divider bank 78. This DFT of the 
system response for the adaptive filter structure 630 is supplied to the similar-DFT-terms 
multiplier bank 632 to be multiplied term-by-corresponding term by the DFT of delayed 
adaptive filter structure 630 input signal as supplied from the analyzer filter 631. The 
synthesizer filter 633 responds to the product output signal from the similar-DFT-terms 
multiplier bank 632 to generate the substantially-echo-free, equalized output signal of the 
adaptive filter structure 630, which is supplied to the circuitry 64 that performs the digital 
synchrodynes to generate an in-phase baseband signal I and a quadrature-phase baseband 
signal Q. 

FIGURE 21 shows a modification of the FIGURE 9 DTV signal receiver to 
replace the transversal adaptive filter 63 by an alternative adaptive filter structure 630. 
The modifications of the FIGURE 9 DTV signal receiver in FIGURE 21 are similar to 
the modifications of the FIGURE 10 DTV signal receiver in FIGURE 20. 

As noted above, when adaptive FIR filtering is used for echo suppression and 
channel equalization, a problem arises because signals with echoes from times previous 
or subsequent to the current time are weighted and combined with the current signal in 
order to suppress the echoes attending it, rather than echo-free signals being so weighted 
and combined. As noted above, further convolution procedures carried out in the DFT 
domain can be used to generate corrections in the weighting coefficients computed for 
the kernel of the adaptive FIR filtering, so performance is improved especially when 
multipath distortion is severe. An analogous problem arises in regard to the weighting 
coefficients supplied to the bank 118 of digital Iowpass filters in FIGURE 20 and in 
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FIGURE 21. The solution to this analogous problem is provided by further convolution 
procedures carried out in the DFT domain, which further procedures are analogous to 
those used for correcting the weighting coefficients computed for the kernel of the 
adaptive FIR filtering. 

5 The equalization methods described in this application can be assisted by 

monitoring the received signal to determine when known data patterns occur and then 
using better estimates of the transmitted signal based on those known data patterns when 
they occur, rather than using estimates of the transmitted signal derived from the response 
of the filtering used for equalization and echo-suppression. By way of example, known 

10 data patterns occur in the initial data segment of each data field in VSB DTV signals 
transmitted in accordance with the ATSC Digital Television Broadcasting Standard 
published in 1995. 

The lowpass filtering of the weighting coefficients of the adaptive transversal 
filtering, which is done to suppress the effects of noise accompanying the received signal, 

15 can be adjusted dependent on how serious that noise is likely to be. If the adaptive 

filtering uses lowpass filtering of the DFT of its system response instead to suppress the 
effects of noise accompanying the received signal, this lowpass filtering can also be 
adjusted dependent on how serious that noise is likely to be. Measurement of the energy 
resulting from differential comparison of the output signal from the adaptive filtering 

20 with the response of a data slicer to that output signal can be used for reckoning the 

strength of noise accompanying a received VSB DTV signal. If there is very little energy 
resulting from the differential comparison, there is not very much noise in the received 
VSB DTV signal and the time constant of the lowpass filtering is preferably shorter in 
nature. If there is considerable energy resulting from the differential comparison, 

25 supposing that the reception channel has been tuned to for some time, there is probably 
substantial noise in the received VSB DTV signal and the time constant of the lowpass 
filtering is preferably longer in nature. 

The DFT methods for adapting the filtering used for equalization and 
echo-suppression that are disclosed earlier in this specification all compare the received 
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signal with estimates of the transmitted signal derived from the response of that filtering 
to the received signal. The estimates of the transmitted signal are accordingly delayed by 
the latent delay of the filter. An alternative way of estimating the transmitted signal is to 
assume that over a long time its frequency spectrum has the ideal channel response. In 
5 the short term this is not generally the case, but if the weighting coefficients are each 
lowpass filtered over a sufficiently long time they should approach desired value. A 
problem with this alternative approach is that the latency of the lowpass filtering of the 
weighting coefficients introduces lag in adapting the filter kernel, which lag is usually 
longer than the latency of the adaptive filter. 

10 The structures described above are designed for synchronous operation at the 

sample rate used in the adaptive filter. At least insofar as the suppression of static 
multipath distortion is concerned, the procedures work better if the samples considered in 
the DFT procedures are in a reasonably large time window extending over 8192 symbol 
epochs, for example. Considering such a large time window, extending several times 

15 longer than the longest delay of echoes with significant energy, facilitates accuracy of the 
DFT calculations for characterizing the channel by reducing the significance of three 
sources of error. Considering such a large time window reduces the effects of the 
post-echoes of symbols preceding the window which post-echoes fall within the window; 
it reduces the effects of the pre-echoes of symbols succeeding the window which 

20 pre-echoes fall within the window; and it reduces the effects of the later symbols within 
the window having their post-echoes fall outside the window. Most of the echoes within 
the DFT time window are of samples also within the window. Considering such a large 
number of samples, however, slows the rate at which the adaptive filtering weighting 
coefficients can be initialized or can react to a sudden radical change in multipath 

25 reception. This is because of the latency of the adaptive filtering with regard to sweeping 
out old symbol values from a time before the reception of new signal symbols and 
because of the further time to supply adaptive filtering response to enough new signal 
symbols to fill the DFT time window. Generally, if the DFT window extends over a 
certain number of symbol epochs, initialization or a correction for changing multipath 

30 conditions will be delayed by at least twice the duration of the DFT time window. The 
delay will be even longer if lowpass filtering of adaptive filtering coefficients is done. A 
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DFT time window extending over 8192 symbol epochs will take a little more than 1.5 
milliseconds to acquire the information needed to begin initialization or correction for 
abrupt change in multipath conditions. This is only a fraction of the time auto-regressive 
adaptive filtering methods require, about 25 milliseconds being the fastest reported time 
for converging an auto-regressive adaptive filter used for equalization and 
echo-suppression of 8-VSB DTV signals. Reducing the duration of the DFT time 
window can reduce the time to acquire the information needed to begin initialization or 
correction for abrupt change in multipath conditions, but the accuracy of the DFT 
calculations for characterizing the channel suffers unless something is done for reducing 
the errors in those calculations. 

The reduction of errors in the DFT calculations for characterizing the channel 
when the DFT time window is shortened is more easily implemented using a computer 
that permits asynchronous operation, lifting the constraint of synchronous operation at the 
sample rate used in the adaptive filter. Then, any of various improved methods can be 
pursued. 

In one improved method of operation, a block of 512 symbols estimated to have 
been transmitted for reception immediately preceding the DFT window for actually 
received signal is zero-extended 512 symbol epochs thereafter and is convolved with the 
just previous adaptive filtering kernel to estimate the post-echoes of these symbols falling 
into the DFT window. These estimated post-echoes can be subtracted from the initial 
512 symbols in the DFT window to reduce the effects of the post-echoes of symbols 
preceding the window on the ensuing DFT calculations. Since post-echoes are the 
principal problem in many multipath reception situations, this can achieve acceptable 
results despite shortening the DFT window to 4096 symbol epochs or so. 

In another improved method of operation the DFT window for actually received 
signal, as shortened to 4096 symbol epochs, is partitioned into a 512-symbol-epoch initial 
section, a 3072-symbol-epoch middle section and a 512-symbol-epoch final section. The 
estimates of a number of the transmitted symbols up to the middle section, including all 
those in the initial section, are zero-extended 512 or so symbol epochs thereafter and are 
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convolved with the just previous adaptive filtering kernel. The convolution results are 
subtracted from the 512 symbols in the initial section and the symbols at the beginning of 
the DFT window to suppress the symbols in the initial section of the window, together 
with their echoes in the initial section of the window, and to suppress their post-echoes in 
5 the middle section of the window. This leaves the initial section of the window occupied 
principally by the pre-echoes of the initial symbols in the middle section of the window. 

The estimates of the transmitted symbols in the final section are preceded by a 
zero-extension of 128 or so symbol epochs and are convolved with the just previous 
adaptive filtering kernel. The convolution results are subtracted from the 5 12 symbols in 
10 the final section and the symbols at the conclusion of the DFT window to suppress the 
symbols in the final section of the window, together with their echoes in the final section 
of the window, and to suppress their pre-echoes in the middle section of the window. 
This leaves the final section of the window occupied principally by the post-echoes of the 
final symbols in the middle section of the window. 

1 5 The DFT window for adaptive filtering input signal, as subjected to the 

corrections per the preceding two paragraphs is divided term-by-corresponding-term into 
the DFT of the convolution of the estimates of the transmitted symbols for the middle 
section of the input signal with the ideal channel response. The result of this procedure is 
the DFT of the adaptive filter kernel. 

20 The adaptive filtering procedures thusfar described have multiplieatively applied 

weighting factors to differentially delayed input signal samples in order to determine the 
corrections that have to be added or subtracted to each signal sample. An alternative 
approach is to apply weighting factors, not to the differentially delayed input signal 
samples, but rather (insofar as possible) to corresponding samples that estimate the 

25 actually transmitted signal. This eliminates a source of noise in the adaptive filter 
response, which noise can interfere with symbol decoding. At the time that adaptive 
filtering is updated, there is the possibility that later-occurring ones of these 
corresponding samples that estimate the actually transmitted signal may not be available 
yet. 
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If the adaptive filter is updated every symbol epoch, only the samples that 
estimate the actually transmitted signal generating post-echoes are available for 
weighting. Supposing the adaptive filter to be configured as a transversal filter, applying 
weighting factors to the samples that estimate the actually transmitted signal to suppress 
5 post-echoes will convert a finite-impulse-response (FIR) filter to a cascade connection of 
a smaller FIR filter for suppressing pre-echoes and an infimte-impulse-response (IIR) 
filter. This is a configuration familiar in auto-regressive adaptive filtering. 

Embodiments of the invention that use auto-regression for reducing error in the 
weighting coefficients, rather than simply lowpass filtering the weighting coefficients to 

10 suppress noise, have been simulated with considerable success. These simulations are 
described in the appendix titled "DESCRIPTION OF SIMULATION WORK". The 
characterization of the channel and calculation of initial weighting coefficients for the 
adaptive filter using DFT methods avoids the start-up problems of "blind equalization" 
auto-regressive filter adaptation methods. Thereafter, the incremental adjustment of the 

1 5 weighting coefficients using auto-regressive methods facilitates the tracking of dynamic 
multipath variations. 

Thusfar, some very important phasing considerations have not been described in 
detail. In fact there are two phases that one needs worry about. A first concern is the 
"demodulator" phase of the synchrodyning or other synchronous detection procedure 

20 used for demodulating the VSB DTV signal to baseband. A second concern is the phase 
of the analog-to-digital conversion as referred to baud phase — that is, to the Nyquist rate 
with respect to the peak of the "desired" principal reception signal at baseband — which 
is sometimes referred to as "symbol clocking phase". If this is off, the data slicer phasing 
is incorrect and the sin(x)/x property of the ideal channel does not have its zeros at the 

25 right spots, even if there are no echoes. This undesirably generates pseudo-echoes in the 
real-only baseband signal. 

If there is the possibility of large echoes being received, there is no assurance of 
accurate measurement of either type of phase until the adaptive filtering substantially 
completes its equalization and de-echoing process. The phase response of the adaptive 
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filtering varies appreciably during equalization and de-echoing process; and, if there is 
appreciable multipath distortion in the actually received signal, the final phase response 
cannot be reckoned directly from the actually received signal as supplied to the adaptive 
filtering. 

5 If the complex adaptive filtering is used, the accuracy of the demodulator phase is 

of little concern so long as its phase rotation with respect to the "desired" principal 
reception signal is not varying at a rate higher than the adaptive filtering can follow. The 
complex adaptive filtering simply corrects for these demodulator phase errors. The 
symbol clocking phase can be simply maintained by a bright-spectral-iine method. For 

10 example, the real and imaginary components of the adaptive filtering response are 
squared and summed to detect the square of the envelope of the "desired" principal 
reception signal The symbol frequency component of this squared envelope is extracted 
by a narrowband FIR filter and is used for generating an automatic frequency and phase 
control (AFPC) signal for the master oscillator determining the timing of sampling in the 

1 5 synchronously operated digital circuitry, including the clocking of analog-to-digital 
conversion. 

Full complex discrete Fourier transforms and inverse discrete Fourier transforms 
are employed in complex adaptive filtering for equalization and echo suppression. If 
real-only adaptive filtering is used, real-only forward and assumed-Hermitian inverse 
20 transforms are used, as done in the simulations described in the appendix. 

If real-only adaptive filtering is used, it can compensate for demodulator phase 
error ranging between -30° and +30°, and it can compensate for errors in symbol clock 
phasing as large as 1/4 to 1/3 of a symbol epoch, lag or lead. A well-designed tracker can 
keep demodulator phase error within the -30° to +30° range as long as dynamic 
25 multipath distortion changes at a reasonably slow rate, the limitation of the tracker in 

tracking dynamic multipath distortion probably not being as stringent as the limitation of 
the adaptive filtering in tracking such distortion. If real-only adaptive filtering is used 
that samples at baud rate, symbol clocking phase can be determined from the responses 
of variously phased matched filters for the PN511 training sequence in the initial data 
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segment of each data field. This requires the master oscillator determining the timing of 
sampling in the synchronously operated digital circuitry to have a stability of better than 
one part per million, which can be obtained using a crystal-controlled oscillator. 

If real-only adaptive filtering is used that samples at twice baud rate rather than 
5 baud rate, the tracker for reducing demodulator phase error can more closely track rapidly 
changing dynamic multipath distortion. Furthermore, symbol clocking phase can be 
adjusted responsive to AFPC signal generated from the data segment synchronizing 
signal or can be continuously adjusted by the Qureshi method. The stability requirement 
on the master oscillator determining the timing of sampling in the synchronously 
1 0 operated digital circuitry is much reduced from one part per million. 

While adaptive filtering used for equalization and echo-suppression in DT V 
signal receivers has been specifically described in this specification and its accompanying 
drawing, the adaptive filtering technique has far-reaching application, as will be 
appreciated by one competent in the adaptive filtering art and acquainted with 
1 5 specification and its accompanying drawing. This application reaches beyond digital 
communication receivers and to adaptive filtering in general. 

DESCRIPTION OF SIMULATION WORK 

20 The simulations done before 30 March 2000 used 4-VSB signal, rather than the 

8-VSB signal used for DTV broadcasting. The use of 4-VSB signal avoided the 
simulations having to take into account the complication of trellis coding. The 
simulations particularly described herein used synchronous equalization with sample rate 
through the filter being baud rate or Nyquist rate, although higher-sample-rate 

25 simulations have also been performed. In the following description of the simulation 
work, M is the number of pre-echo terms being considered, and N is the number of 
post-echo terms being considered. The word "data" refers to digital baseband signal as 
recovered by demodulating 4-VSB (or 8-VSB) DTV signal, not to the results of symbol 
decoding the digital baseband signal. 
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Receivers that use FIR filters for suppression of pre-echoes must compute, 
iteratively or otherwise, special filter coefficients which differ from the measured echo 
ensemble. The DFT methods used in the simulations suppressed echoes using the echo 
ensemble itself, directly as measured, in filtering procedures more similar to the IJR 
filtering used in the prior art for canceling post-echoes. In these simulations, however, 
these filtering procedures were used to suppress pre-echoes as well as post-echoes. An 
echo-measurement register was provided for temporarily storing the measured echo 
ensemble, as initially determined during start up, and as subsequently updated by 
echo-tracking procedures. This echo-measurement register was separate from an 
echo-suppression register used for accumulating the results of echo-suppression 
procedures. This echo-measurement register contained enough actual storage locations to 
store the range of echos to be corrected and contained additional storage locations to store 
zero-extension information. These additional storage locations could be virtual storage 
locations, rather than actual ones. 

The simulations done before 30 March 2000 employed discrete Fourier 
transforms (DFTs) with 4096 terms for echo-suppression computations. DFTs of blocks 
of data points, which blocks were 3072 = 4096 * (3/4) symbol epochs long, provided the 
information used to remove pre-echoes. Post-echoes were removed using the previously 
sliced data. The echo-suppression register had 4096 storage locations, each associated 
with a respective one of successive symbol epochs of the signal to be equalized. As an 
initial step in the processing of each block of input data points, a guard interval of 5 12 
zeroes was inserted into the final section of an array of 4096 storage locations within the 
echo-suppression register corresponding to symbol epochs occurring after the block of 
data points that was being processed at the time. The block of incoming data 3072 
symbol epochs long was loaded into corresponding ones of the 3072 storage locations in 
the middle section of the echo-suppression register. The 512 storage locations in the 
initial section of the echo-suppression register, containing storage locations 
corresponding to symbol epochs occurring before the block of data points that was being 
processed at the time, were initially filled with the sliced data of the corresponding 
elements from the previous block. 
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Then, using the measured echo ensemble read from the echo-measurement 
register, the echoes of the first point into the next 51 1 points were calculated; the echoes 
of the second point into the next 510 ones were calculated, the echoes of the third point 
into the next 509 points were calculated, etc. This broke up the propagation of echoes 
5 from previous data so that the short 512 zero pad at the end sufficed. This procedure 
required (M+N)*(M+N+l)/2 multiply-and-add steps, where was the number of pre-echo 
terms and N the number of post-echoes. However, these calculations had already been 
done and stored up during the processing of the previous block. The discrete Fourier 
transform of the input data points in the current block was then computed. 

10 The echo-measurement register had 4096 apparent storage locations for 

temporarily storing the echo ensemble, as measured over a period of 512 symbol epochs 
or so, and zero-extended to fit into a 4096 symbol-epoch-long block. The discrete 
Fourier transform of this zero-extended echo measurement was calculated. Then, the 
DFT of the input data points was divided term-by-corresponding-term by the DFT of the 

1 5 echo measurement. The inverse transform of the resulting DFT was the directly 
de-convolved data, which could have been simply sliced and used as they were. 

Instead, in order to reduce the effects of noise, an array 4096 symbol epochs long 
was generated, with the last 4096 - 512 = 3584 elements set to zero and with the first 512 
elements set to the terms of the echo which represented the pre-echo only. The values for 

20 the "desired" term and all the post-echoes were set to zero. This array was then discrete 
Fourier transformed. The pre-echo-only transform was then multiplied, 
term-by-corresponding-term, by result of the division process described in the previous 
paragraph. The resulting DFT was then inverse-transformed. The final result was an 
array containing 3072 values presumed to be useful, which values provided an 

25 approximation (containing noise and errors due to mis-measurement of the echo 
ensemble) to a signal that solely consisted of the pre-echoes. There was no principal 
signal in this final result, nor any post-echoes. 

The operation continued as follows. The value of each data point to be sliced was 
calculated by taking each incoming data point, first to last in the 3072 long block stored 
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in the echo-suppression register, and subtracting from it the corresponding pre-echo sum 
created in the DFT process described above, thereby updating the corresponding 
locations in the echo-suppression register. Then, in a further updating of the 
corresponding locations in the echo-suppression register, each of the post-echo terms was 
removed by subtracting, one by one, the product of the echo for that distance in the past 
and the already sliced value of the data point for that time. This required N 
multiply-and-subtract operations, where N was the number of post-echo terms used. 
Each multiplicand, being supplied from the data slicer, was only three bits. The 
foregoing procedures generated the de-echoed data point stored at a respective storage 
location in the echo-suppression register, which de-echoed data point was then sliced or 
actually run through the Viterbi decoder. This procedure was repeated until each and 
every one of the 3072 data points in corresponding storage locations of the 
echo-suppression register was used to update the echo-suppression register for removing 
any pre-echoes of that data point and then to further update the echo-suppression register 
for removing any post-echoes of that data point. 

The foregoing procedures were used during initialization or after a radical change 
in multipath reception conditions. During normal multipath reception conditions, these 
procedures were augmented by a feedback mechanism that reduces error in the 
echo-supression calculations. This feedback mechanism proceeded together with the 
process of removing the echoes from the data points and slicing them, but the operation 
of this feedback mechanism was delayed in regard to when each data point is sliced. A 
delay of M points provided time for the pre-echoes of any given "desired" point to be 
calculated, after which the feedback mechanism operated with respect to each data point 
as follows. The error term that was subtracted from this data point M points earlier to 
remove the pre-echoes was temporarily stored and then added back in to recover the data 
point with the original pre-echoes in it, updating the corresponding storage location in the 
echo-suppression register. Then, in M multiply-and-subtract operations using the 
corresponding echo points and the M sliced values generated since this data point was 
sliced, the pre-echo was subtracted using the sliced data point instead of the DFT 
estimate to update the corresponding storage location in the echo-suppression register. 
This produced less noise and instability in a new estimate of the de-echoed data point 
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stored at that corresponding storage location in the echo-suppression register. The 
corresponding sliced (or quantized) data point was subtracted from this new estimate, 
generating the error in the calculation of the data point. Then, in the simplest form of the 
simulation, the product of this error times each of the M+N+l sliced data points 
corresponding to N post-echoes, the M pre-echoes and the main signal was calculated. 
The resulting products were each multiplied by a constant attenuation factor u and added, 
one by one, to the ensemble of echo points temporarily stored in the echo-measurement 
register. This procedure calculated the correlation function of the error and the sliced 
data, then applied it as the correction, described above. This procedure used the full 
available information to update the echo. This procedure entailed many multiply-and add 
operations, (M+N+l) of them, but note that they were not the full precision of the 
original data. Each sliced data point had only 3 bits; and the error had a maximum value 
of+1 or-1 (assuming sliced levels of -7, -5, ... 7), because any value greater than that is 
likely to be an error. Note that in the process of removing the pre-echo using the sliced 
data points there was also generated, as an intermediate result, the exact terms needed in 
the process of generating the data described above for use as input to the second DFT for 
the next block. 

This is the basic description of the simulations done before 30 March 2000. Care 
was taken not to allow the quotient of the two DFTs to go out of range because of the 
divisor being too small. Note that nothing was done to the transforms while in the 
frequency domain, except operating on each element exactly the same way. This meant 
that there was no need to do the normal Cooley-Tukey "shuffle" of the elements to get 
them into frequency order for the forward transform and out of that order for the inverse 
transform. Normally, these DFTs would be computed in dedicated hardware hard-wired 
to take advantage of the symmetries and known zero elements involved. 

The importance of the weighting coefficient u. was noticed when performing the 
simulations. The larger u. is, the more rapidly that dynamic multipath can be followed. 
The smaller u- is, the more stable the process of following the echo ensemble is. 
Preferably, the noise in the echo and the rate of speed of its change are monitored and the 
value of u is adjusted accordingly. In general, values of ix from 1/800 to 1/160,000 are 
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typical. The digital AGC is another coefficient in the system that should also be updated 
continuously. 

There are two limitations on the tracking of dynamic multipath using this method. 
A first limitation is that the value of jx cannot be too large, or the errors in the echo 

5 measurement will cause a blowup. However, if jx is too small, the echoes cannot be 
followed fast enough. This can to some extent be helped by using a more complicated 
method for updating the echo than the simple addition described above. What is 
described produces a simple exponential weighting of previous values of the echo, falling 
to 1/e in 1/jx samples. More complicated differential equation solvers instead of the 

10 simple addition will change the properties of this. A second limitation on the tracking of 
dynamic multipath is the fact that, if there are pre-echoes, then part of the de-echoing of 
them applied before slicing (but not to the feedback) is done stepwise every 285 
microseconds. This of course results in these values being in error somewhat part of the 
time. 

15 Provisional U. S. patent application serial No. 60/193,301 indicates that one can 

use smaller than 3072 long blocks for intermediate tradeoffs between hardware 
requirements and robust dynamic echo tracking. Simulations were performed using 3072 
long blocks in order to emulate the dynamic multipath performance of 2K COFDM, 
Such simulations would conveniently provide roughly equivalent performance to 2K 

20 COFDM for pre-echoes (and, of course, better performance for post-echoes.) Besides, 
when the simulations before 30 March 2000 were done, processing 3072 data points at a 
time in the DFT procedures was viewed as being convenient. A block of 3072 data 
points was considered not to be too small, so the necessary padding to prevent 
wraparound did not become too large a fraction of the DFT. It was considered not to be 

25 too large, so that the echo tracking would be poor because the tracking of pre-echoes 
being harmed. 

Since 30 March 2000, further simulations have revealed that processing 2048 
symbol epochs at a time is a more practical upper limit for echo tracking procedures 
using 4096-point DFTs. Processing blocks of more than 2048 symbol epochs results in 
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loss of tracking capability when Rayleigh echo ensembles are received Processing of 
fewer data points at a time permits better tracking of fast-changing dynamic echoes. 
Echo tracking procedures using 2048-point DFTs should not process more than 256 
symbol epochs at a time, in order to avoid loss of tracking capability when Rayleigh echo 
5 ensembles are received. 

The method used for starting up the foregoing simulation procedures depended on 
using the PN5 1 1 training sequence. The analog portion of the DT V receiver was 
presumed to have achieved lock. While waiting for the training sequence to arrive, the 
DFT of the autocorrelation function of the incoming random data was computed using 

10 either the FIR filters described above or the block DFT hardware also used above. Using 
a total of perhaps 50,000 data points was considered to be optimal. Then, the PN51 1 
training sequence was searched for, by calculating correlation functions between the 
incoming data and a stored copy of the PN5 1 1 training sequence, using DFTs. Once the 
PN51 1 training sequence arrived, this cross-correlation was available as a measurement 

15 of the echo and could be used to start the simulation procedures described above. 
However, the measurement of the echo afforded by this cross-correlation was 
considerably corrupted by noise, part of which noise arose from post-echoes of the 
random data that preceded the PN 51 1 sequence. These post-echoes superposed 
themselves on the PN 51 1 sequence they overlapped in time. 

20 A preferable procedure was to use the cross-correlation of the known PN5 1 1 

training sequence and the data to align the known PN5 1 1 training sequence and the data, 
after which the known PN5 1 1 training sequence as properly scaled and filtered was 
subtracted from the data. The filtering of the PN5 1 1 training sequence included 
convolution with the ensemble of strongest echoes as determined from the 

25 cross-correlation. Presumably, this could be done using the hardware used in the running 
process described above. Subtracting echoes as small as 0.08 of the "principal" signal 
can be useful, but even doing only down to 0.2 was found to be very helpful. Then, using 
DFTs, the cross-correlation of the "reduced" data and the 511 point sequence was 
computed. This cross-correlation had the principal peak and the principal echoes 

30 previously subtracted missing, so they were added back in. However, the resulting sum 
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had a much smaller error than the original cross-correlation of the known PN5 1 1 training 
sequence and the data. This is because a 51 1 -element pseudorandom noise (PN) 
sequence is simply too short to be a truly complete "pseudorandom" sequence and the 
cross-correlation of the sequence with itself has a large residual. The DFT of that 
5 cross-correlation was then computed. At this point, the real part of the DFT of this 
cross-correlation was replaced by the real part of the much longer averaged 
autocorrelation of the random data calculated earlier, which had a smaller average error. 

This was possible because of the fact that the autocorrelation of the data has the 
correct amplitude spectrum, but an incorrect phase spectrum. This modified DFT was 
1 0 then inverse-transformed to give the echo. This last step reduced the error typically by 
about the square root of two, while the iteration procedure on the cross-correlation had a 
much larger effect. This start-up procedure took some time, so the incoming data 
occurring during the time it took was stored up; and then, in order to catch up, the clock 
for the continuing de-multipath procedure was run faster than otherwise needed. 

1 5 Finally, there were two empirical procedures that were used for startup and 

recovery from "bad channel" conditions that did not quite cause loss of lock. In the first 
of these empirical procedures, a running measure of a histogram of the error in each 
"eye" was calculated from data used in the feedback loop. Then, every hundred points or 
so, there are measurements made of the depths of the valleys between the peaks (eight 

20 peaks in 8-VSB, but only four peaks in the 4-VSB used in the simulations made before 
30 March 2000). If the valleys were larger than one-sixth the peaks, the feedback system 
was modified in the following way. If the error was more than one half the distance 
between two slice positions (i.e., one for 4-VSB with levels -3, -1, 1, 3), then it was not 
added into the feedback. This prevented wrong slices from causing feedback in the 

25 wrong direction. This was modified in later simulations using 8-VSB signals, where 
prediction based on the state of the Viterbi decoder was used to do a similar thing. 

The second empirical procedure was only used in the start-up procedure. If one 
trial failed, on the next training sequence that came along the measured echo was 
truncated, so that terms with absolute value no larger than 0.04 of the biggest peak were 
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set to zero. This procedure was unstable. So, as soon as the eyes opened appreciably, the 
procedure was discontinued. 

The feedback mechanism used in the simulations required (M+N) 
multiply-and-add operations per data point. This can be reduced dramatically with a very 
substantial reduction in dynamic multipath performance by the following alternative 
method. The error signal, measured as described above, is collected up for all 3072 
points until the end of the block. The result is zero-extended to 4096 samples and 
discrete Fourier transformed. The 3072 sliced data points are zero-extended to 4096 
samples and discrete Fourier transformed. The former DFT is either divided by the latter 
DFT, or is multiplied by the complex conjugate of the latter DFT and inverse transformed 
to generate the correction to be applied to the echo. Note that there are two possible ways 
to perform this alternative method, one using a de-convolution step and the other using a 
correlation step, and they seem to have equal performance. This alternative method 
requires less hardware to perform than the one previously described, but has to be done 
after the block is finished, increasing the speed at which the processor must operate. It 
also has the disadvantage that it provides much less frequent updates of the echo 
ensemble measurement, which updates are applied in block form. This means it will not 
follow echoes as fast. As of 30 March 2000, it appeared that the simple addition of a 
prescribed fraction of the measured correction does not perform adequately. Instead, a 
more complicated differential equation solver type system appeared to be required, 
perhaps one using non-linear error change prediction methods based on stored up 
previous echo ensembles. 

The simulation methods thusfar described provided full tracking for post-echoes, 
but not for pre-echoes. Provisional U. S. patent application serial No. 60/193,301 
proposed a simulation method that modifies the above-described echo-suppression 
method to treat pre-echos much faster, which modified method should have identical 
results for echo ensembles consisting only of post-echoes. The proposed method uses a 
feedback mechanism identical to that in the first-described simulation method and works 
as follows. The proposed method nominally works on blocks of data that are the post- 
echo length, N symbol epochs, long. At the start of processing each block, as before, 
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data points have been sliced (actually Viterbi decoded) up to the start of the block. The 
echo has also been tracked to this point. Discrete Fourier transforms at least 2N long are 
employed to avoid wraparound. The first N elements of an array are filled with the last N 
known sliced data points (array A). Then the end of the array is padded with zeros as 

5 data points to avoid wraparound. The post-echo is placed into an equal-length array B 
padded with zeros at both ends (see later). The arrays A and B are Fourier transformed to 
generate arrays C and D, respectively. The arrays C and D are multiplied together, 
element-by-corresponding-element, and the result is discrete Fourier transformed to 
generate an array E. The array E comprises the echoes of the last N sliced points 

10 propagated into the future up to N. The next N incoming data points are placed into an 
array of the same length, starting at the position next in that array to where the last "old" 
sliced data point was in the array A and the forward-propagated post-echoes in array E 
are subtracted away to form an array F. This is then transformed to generate an array G. 
The pre-echoes are then placed into a same-size array H so that they occupy the correct 

1 5 relative position in the array H relative to where the post-echoes were before in the array 
B, which is the reason the post-echoes needed padding at the beginning of the array B. 
The array H is discrete Fourier transformed to generate an array J. The array G is 
divided, element-by-corresponding-element, by the sum of the arrays D and J. The 
resulting array K is the DFT of the completely de-convolved data of the N points. The 

20 array K is multiplied, element-by-corresponding-element, by the array J to generate an 
array L consisting of the DFT of the pre-echoes only. The array L is inverse- 
transformed to generate an array M. 

Now, the method proceeds using the array E, the array M, and the original 
untreated data for this block. The array M is subtracted from the incoming data resulting 

25 in an array N with the pre-echoes removed. Array E is subtracted, 

element-by-corresponding-element, from the array N to generate an array O with the first 
element being completely de-echoed, the second needing to have the echo of the first 
removed, the third needing to have the echoes of the first and second removed, etc. This 
process of removal of these post-echoes is done an element at a time using recently sliced 

30 (actually Viterbi decoded) numbers and the current echo measurement and subtracting 
element by element This last part requires N*(N+l)/2 additions and subtractions. After 
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each point is in turn rid of its post-echoes, it is sliced and used as input to the next step. 
This is similar to the first method. The difference is that here a part of the removal of 
post-echoes from the current array is done by transform methods. 

An important point about this method is that it uses much smaller blocks, so that 
5 one can implement feedback to an analog phase-locked loop much more easily. One 
would presumably need blocks 1024 long. A second point is that it uses much more of 
the information from sliced data than the original method to remove pre-echoes and so 
should be more reliable. Blocks longer than N could be used, in order to waste less of the 
1024 long transforms, but this would lose some of this advantage. The proposed method, 
1 0 in terms of quality of pre-echo removal, would look more and more like the method used 
in the original simulations, and would eventually be the same except for the few elements 
at the beginning of the block. 

A careful examination of this proposed method reveals that it needs a total of .s ix: 
transforms based on a naive implementation. And the four forward transforms all have 

1 5 more than half of the input zeros. This appears to be inefficient, but this is not the case. 
There are techniques that allow doing both halves at once and then sorting out the results 
at the end. So, the arrays C and D are generated concurrently. So are the arrays G and J. 
See as a reference, for example, "Numerical Recipes in C" (or Fortran). One would not 
need, also, to do the "shuffle" normally done if one used a decimation in frequency 

20 algorithm. This is because there are no operations dependent on the actual frequency, 
just multiplication and divisions on like elements. 

The simulations suggest that the strongest reception-path response is not always 
the best choice for the "desired" principal reception signal. The simulations indicate that 
it is better to have a 1 50% post-echo than a 67% pre-echo. If two of the strongest 
25 reception-path responses are reasonably close in time and within a factor of 1 . 5 in peak 
amplitude, it is simpler to chose the earlier reception-path response as the principal 
reception signal and the later reception-path response as a post-echo, then let the adaptive 
filtering proceed. It will work fine. If that is not true, since a tracker is used that can 
catch up at least a bit if it gets behind on occasion, this is one of those occasions. The 
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control computer just decides anew which is the new "principal peak" and directs 
reloading of all the registers, backward in time if necessary, and directs that start-up 
procedures be employed again. This is a messy procedure, but is a problem that VSB has 
in common with COFDM. 
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